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ABSTRACT 

In the search for life on Earth-like planets around other stars, the first (and likely only) 
information will come from the spectroscopic characterization of the planet’s atmosphere. Of 
the countless number of chemical species terrestrial life produces, only a few have the distinct 
spectral features and the necessary atmospheric abundance to be detectable. The easiest of these 
species to observe in Earth’s atmosphere is O2 (and its photochemical byproduct, O3). But O2 
can also be produced abiotically by photolysis of CO2, followed by recombination of O atoms with 
each other. CO is produced in stoichiometric proportions. Whether O2 and CO can accumulate 
to appreciable concentrations depends on the ratio of far-UV to near-UV radiation coming from 
the planet’s parent star and on what happens to these gases when they dissolve in a planet’s 
oceans. Using a one-dimensional photochemical model, we demonstrate that O2 derived from 
CO2 photolysis should not accumulate to measurable concentrations on planets around E- and 
G-type stars. K-star, and especially M-star planets, however, may build up O2 because of the low 
near-UV flux from their parent stars, in agreement with some previous studies. On such planets, 
a ‘false positive’ for life is possible if recombination of dissolved CO and O2 in the oceans is slow 
and if other O2 sinks (e.g., reduced volcanic gases or dissolved ferrous iron) are small. O3, on 
the other hand, could be detectable at UV wavelengths (A < 300 nm) for a much broader range 
of boundary conditions and stellar types. 

Subject headings: planets and satellites: atmospheres - planets and satellites: terrestrial planets - planet- 
star interactions - ultraviolet: planetary systems 

- 1. Introduction 
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ing ground-based radial velocity measurements 
(Anglada-Escude et al. 2012 1 and transit obser¬ 
vations from NASA’s Kepler Space Telescope 


Borucki et al. 

2012 

2013 

Quintana et al. 

2014 


has led to increased interest in the question 
of whether other Earth-like planets exist and 
whether any of them might be inhabited. If 
Earth- and super-Earth-sized planets can be found 
around nearby stars with the Transiting Exoplanet 
Survey Satellite (TESS) (Ricker et al. 20141, we 
may be able to analyze their atmospheres with 
transit spectroscopy using the upcoming James 


Webb Space Telescope (JWST) (Deming et al. 


2009 1 . The next generation of space telescope 


mission concepts, such as the Advanced Technol¬ 
ogy Large Aperture Space Telescope (ATLAST) or 
the Large Ultraviolet/Optical/Infrared Telescope 
(LUVOIR) (e.g., [Postman et al.|2010[ France et al. 


20151, could use reflection spectroscopy to char¬ 


acterize many more such planets. Future large 
(30-40 m aperture) ground-based telescopes could 
also characterize some of the nearest detected ter¬ 
restrial planets (Kawahara et al. 2012[ Snellen 


et al.|[20T3|. More focused mission concepts, such 


as NASA’s Terrestrial Planet Finder (TPF-C) or 
Webster Cash’s New Worlds Observer Mission 
(also known as TPF-0), might also characterize 
a large number of Earth-like planets. Earth-like 
exoplanets could be studied at longer wavelengths 
using emission spectroscopy by interferometers 
such as NASA’s TPF-I or ESA’s Darwin mission 
(Cockell et al. 2009). None of these missions is 
being actively pursued by NASA at the moment 
but they could be selected for development several 
years from now if they receive a favorable rating 
in the next Astronomy and Astrophysics Decadal 
Survey. Smaller probe-class concepts featuring 
a coronagraph (Exo-C) or a starshade (Exo-S) 
would be a low-cost stepping stone from JWST to 
a dedicated TPF mission, and could fly in concert 
with the next UVOIR space telescope. 


1.1. Disequilibrium pairs of gases as a 
biosignature 

Because atmospheric characterization of exo¬ 
planets can be done simultaneously with direct 
detection, the question immediately arises as to 
what gases, or combination of gases, might pro¬ 
vide evidence for the existence of life on such 
a planet. A capable review of this subject has 


been given by Seager & Bains (2015). We pro¬ 


vide our own brief synopsis here. The literature 
on remote life detection extends back well before 


the exoplanet era. Lovelock (1965) first suggested 


that the presence of oxygen in a planet’s atmo¬ 
sphere, alongside hydrocarbons, would constitute 
a reliable biosignature. Later, [Lippincott et al.| 
(1967) extended this to specific gas pairs such as 


O 2 and CH 4 or N 2 O. These gases are many or¬ 
ders of magnitude out of thermodynamic equilib¬ 
rium with each other, and all have sources that are 


dominantly biological. Lederberg (1965) had al¬ 


ready proposed that extreme thermodynamic dise¬ 
quilibrium, in general, would be good evidence for 
life. This concept has been explored recently by 


Bains & Seage 

r ( 2012 ) (see also 

Seager et al. 

(2012 

2013a|b ) and 1 

Crissansen-Totton et al. 

(2015 

), but 


no additional, potentially observable, disequilib¬ 
rium redox pair has been described. Moreover, 
photochemical modelers (e.g., Kasting|1990 2014 


Zahnle et al. 2008) have shown that CO-rich atmo¬ 


spheres can arise under a variety of circumstances. 
CO is a high free-energy compound that would 
be well out of thermodynamic equilibrium in al¬ 
most any plausible abiotic terrestrial atmosphere. 
But, for this same reason, it should be an excellent 


source of metabolic energy for microbes (Kharecha 


et al. 2005 Kasting 2014), and so its presence in 


a planet’s atmosphere might actually be consid¬ 
ered an anti-biosignature { Zahnle et al.|2008 ). We 
thus agree with Seager & Bains (2015) that the 


use of thermodynamic disequilibrium as a biosig¬ 
nature cannot be easily generalized; however, the 
more specific life detection criterion proposed by 
[Lovelock] and [Lippincott et al.| is still useful. 

Herein lies a difficulty, however. O 2 is highly 
abundant (21% by volume) in Earth’s atmo¬ 
sphere and should be readily detectable by a first- 
generation TPF-type mission by way of its ‘A’ ab¬ 
sorption band at 760 nm ( Owen[[T980 Des Marais 
et al. 2002). But CH 4 is present at only 1.7 


ppmv and would likely not be detectable by a 
first-generation TPF-type telescope, despite the 
fact that it absorbs in both the visible and near- 
IR. Other, clearly biogenic, gases such as methyl 
chloride (CH 3 CI) and dimethyl sulfide (DMS) are 
also unlikely to be detectable by TPF (Domagal- 


Goldman et al.pOll Seager fc Bains[[26l5 ). This 

leads to the question that is the main topic of this 
paper: Can O 2 by itself be considered a biosig- 
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nature and, if so, under what conditions is this 
interpretation valid? 


1.2. O 2 by itself as a biosignature 


1.2.1. Historical development 

We begin by acknowledging that many previous 
researchers have suggested, directly or indirectly, 
that O 2 is not a good bioindicator. In a series 


of papers written during the 1960’s, Berkner &: 
Marshall] (1964 1965, 1966 19671 calculated that 
atmospheric O 2 could build up to as high as 10“^ 
PAL (times the Present Atmospheric Level) as a 
consequence of photolysis of H 2 O followed by es¬ 
cape of hydrogen to space. This amount of O 2 is 
just on the verge of being detectable by TPF, pri¬ 
marily through its photochemical byproduct, O 3 


(Leger et al. 1993 Segura et al. 20031. And, as 


we discuss further below, this is also the O 2 level 
at which confusion with biologically produced O 2 
begins to become an issue. Several years later, 


Brinkmann (1969) predicted prebiotic O 2 levels as 


high as 0.27 PAL, which would certainly be de¬ 
tectable by currently envisioned technology. Both 
studies, however, were performed before the pro¬ 
cess of hydrogen escape from Earth’s atmosphere 
was fully understood. We know now that the 
rate of hydrogen escape is limited by diffusion 


through the homopause, near 100 km (Hunten 
1973 Walker |1977 ), and this diffusion rate, in 
turn, depends on the total hydrogen mixing in the 
stratosphere: ftot(H) = f(H) -|- 2 f(H 2 ) -I- 2 f(H20) 


-I- 4 f(CH 4 ) -I- ... (see, e.g., Kasting & Catling 


20031. The modern Earth’s stratosphere is very 
dry (f(H 2 0) « 3-5 ppmv), and reduced gases are 
also at or below the ppmv level, so the current 
rate of O 2 production from hydrogen escape is rel¬ 
atively slow. 

More recent claims of high abiotic O 2 levels 
have been made by researchers using improved 
photochemical models and more realistic assump¬ 
tions about hydrogen escape. jCanuto et al. (1982 


19831 predicted high abiotic O 2 levels from CO 2 
photolysis by a highly UV-active young Sun. 
Their assumed UV levels were unrealistically high 
for the early Earth, and their model did not obey 
redox balance (more on this below), so their ap¬ 
parent false positive was probably not correct 


(Kasting et al. 19841. Rosenqvist & Chassefiere 


(19951 predicted high abiotic O 2 concentrations, 


but their model also failed to consider redox bal¬ 
ance ( Kasting|1995 1. Selsis et al. (2002) predicted 
high abiotic O 2 levels in dense CO 2 atmospheres, 
but their model suffered from this same problem 


(Segura et al. 2007). Given the critical impor¬ 


tance of redox balance to abiotic O 2 levels and 
the absence of any consideration of it in numerous 
papers in the literature, we include a discussion of 
this topic below. 

1.2.2. Planets that sujfer extensive water loss 

One class of apparent false positives for life 
involves planets that lose a lot of water. A 
good example is a planet interior to the habit¬ 
able zone that loses its water by way of a runaway 


or moist greenhouse ( 

Kasting 

1988 

Wordsworth & Pierrehumbert 

2013) 


porized Earth ocean would have a surface pressure 
of ~270 bar. Thus, if such a steam atmosphere 
was photolyzed by stellar UV, and if the hydrogen 
escaped to space, the planet could be left with a 
240-bar O 2 atmosphere, having lost just over 10 
percent of the ocean mass. But this type of false 
positive is not too worrisome because we would not 
see H 2 O absorption bands, unless we were unlucky 
enough to observe the planet during the brief pe¬ 
riod in which the ocean had evaporated but not 
yet been completely lost; thus, we would likely de¬ 
duce that the observed O 2 was the product of a 
runaway greenhouse. Exogenous water delivered 
later in the planet’s lifetime should be also be 
short-lived, as the planet should be even more sus¬ 
ceptible to water loss when its parent star is older 
and brighter. A similar process could happen to a 
planet in the middle of the HZ if the background 
pressure of noncondensable gases such as N 2 , CO 2 , 
and Ar was small, because H 2 O would then be 


a major atmospheric constituent (Wordsworth & 


Pierrehumbert 2014). Water loss would continue 


until the background gas (in this case, O 2 ) built up 
to ~0.2 bar. This type of false positive would be 
more difhcult to distinguish spectroscopically, but 
might be eventually identified by looking for O 2 - 
O 2 and N 2 -N 2 features that could constrain the 
atmospheric partial pressures of these two gases 


(Misra et al.||2014a 

Schwieterman et al.||2015 

I- 

A recent variant 

on the runaway greenhouse 

concerns planets orbiting M stars. Luger & Barnes| 

(2015) (see also 

Ramirez & Kaltenegger (2014) 


and Tian & Ida (2015)) have shown that because 
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of the long, highly luminous, pre-main-sequence 
lifetime of M stars, planets orbiting within their 
habitable zones could lose most or all of their 
water early in their histories, with accompanying 
buildup of atmospheric O 2 . This would be true 
for terrestrial planets around M stars with long 
pre-main sequence lifetimes, regardless of their 
(largely unconstrained) initial volatile abundances 
(e.g., [Raymond et al. 2004, 2007b|a 


2009 Ciesla 


et al.||2015 ). Such planets should also lack strong 
H 2 O absorption features, much as Venus does to¬ 
day, and so the O 2 signal could probably be dis¬ 
counted as evidence for life (Venus’ 30 ppmv of 
H 2 O is detectable spectroscopically, but only by 
looking at high spectral resolution through near- 


IR ‘windows’ on the planet’s nightside (Meadows 


& Crisp 1996) - Venus’ H 2 O is not detectable 
at TPF-type spectral resolution ( Meadows|2005 1 ). 
Hamano et al. (2013) have suggested that plan¬ 


ets that start off with a runaway greenhouse at¬ 
mosphere could have their O 2 absorbed by the 
magma ocean, perhaps removing the possibility 
of a false positive altogether. But one can imag¬ 
ine M-star planets that are resupplied with water 
at some later time by some process akin to the 
Late Heavy Bombardment in our own Solar Sys¬ 
tem, or that migrate into the HZ after their host 
star evolves onto the main sequence. Hence, we 
should still be concerned about false positive on 
M-star planets, as we discuss further below. 

1.2.3. Planets with frozen surfaces 

A slightly less obvious false positive is a planet 
just beyond the outer edge of the habitable zone, 
like early Mars, but with roughly twice Mars’ mass 
(Kasting 1997 2010). The frozen surface would 


not take up much oxygen once the surface layer 
itself was oxidized, and the planet should be large 
enough to hold onto its atmosphere, but perhaps 
not large enough to maintain volcanic outgassing 
of reduced species like H 2 . In that case, pho¬ 
todissociation of H 2 O in the planet’s atmosphere, 
followed by escape of hydrogen to space, could 
conceivably cause substantial amounts of O 2 to 
accumulate, just as on the Venus-like planet de¬ 
scribed above. The martian atmosphere itself con¬ 
tains ^0.1% O 2 , produced by this process. Indeed, 
Mars’ O 2 concentration would probably be even 
higher except that the planet’s low gravity allows 
oxygen to escape by way of nonthermal processes 


such as dissociative recombination of ions, e.g., 
02 ^-l-e—^O-I-Ojin which the resultant O atoms 
have sufficient kinetic energy to escape. Most of 
these nonthermal loss processes for oxygen would 
cease to operate on a planet twice as massive as 
Mars, allowing O 2 to accumulate almost indefi¬ 
nitely. The planet’s upper atmosphere would be 
dry, making it unlikely to be distinguishable using 
transit transmission spectroscopy, but this type of 
false positive could probably be identified with di¬ 
rect imaging because the planet surface would be 
cold, and hence its lower atmosphere should be 
deficient in H 2 O. 


1 . 2 . 4 . More troubling false positives: low-outgassing 
planets and planets around M stars 

Recently, several researchers have suggested the 
existence of more troubling false positives that 
would be harder to identify. It is these studies that 


motivate the present paper. Hu et al. (2012) con¬ 


cluded that C02-dominated worlds around Sun¬ 
like stars may build up oxygen if the outgassing 
rates of reduced gases (such as methane or hy¬ 
drogen) are small. Their model included an at¬ 
mospheric hydrogen budget, along with realis¬ 
tic assumptions about the rate of hydrogen es¬ 


cape. More recently, Tian et al. (2014) have sug¬ 


gested that C02-rich planets orbiting M stars may 
build up significant amounts of O 2 because of the 
higher ratio of far ultraviolet to near ultraviolet 
(FUV/NUV) stellar radiation. Far-UV radiation 
(A <200 nm) can photolyze CO 2 , producing CO 
and O atoms: 


CO 2 +hiy CO -f O 

Direct recombination of CO with O is slow, how¬ 
ever, because that reaction is spin-forbidden. 
Hence, O atoms recombine with each other to 
form O 2 : 


O-bO-bM ^ 02-bM 

Here, ‘M’ is a third molecule, needed to carry of 
the excess energy from the collision. Recombi¬ 
nation of O or O 2 with CO can be catalyzed by 
the byproducts of H 2 O photolysis, which occurs at 
wavelengths < 240 nm. (All of these photochem¬ 
ical models, including the independent one by Hu 
et ah, extrapolate the measured H 2 O cross sec¬ 
tions from 196 nm out to 240 nm, which is the en- 
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ergy threshold for breaking the H 2 O bond (Kast- 


ing & Walker 19811. The H 2 O absorption cross 
section is low, but finite, at these wavelengths, 
giving rise to substantial H 2 O photolysis all the 
way down to the planet’s surface.) For example, 
one such catalytic cycle (originally proposed by 
Donahue ( |1969 ) and McElroy & Hunten (1970)) 


begins with 


H 2 O + hi/ ^ OH + H 


followed by 


CO + OH 
O 2 + H + M 
_ HO2 + O 

net: CO + O 


CO 2 + H 
HO 2 + M 
O 2 + OH 
CO 2 


The wavelength regions for photolysis of H 2 O and 
CO 2 overlap with each other (Fig. [^; however, 
only H 2 O can be photolyzed longward of 200 nm, 
and so the efficiency of such catalytic cycles de¬ 
pends critically on the flux of stellar near-UV ra¬ 
diation. We have summarized these reactions in 
Fig.d 

Domagal-Goldman et al. suggested that total 
far-UV (<200 nm) fluxes in excess of solar levels 
are responsible for the generation of abiotic O 2 
and O 3 . In contrast, Tian et al. suggest that 
the higher ratio of far-UV (<200 nm) to mid- and 
near-UV (200-400 nm, hereafter referred to sim¬ 
ply as the near-UV) radiation coming from an M 
star (as compared to the Sun) is responsible for 
increased O 2 and O 3 production. For the spec¬ 
tra used in this study, the difference is centered 
at ^170 nm, characterized by a marked change in 
flux for the F, G, and K stars, with little to no 
change in flux for the M stars (Fig. [^. Many M 
stars are highly active and so emit an excess of 
far-UV radiation from their chromospheres. Near- 
UV radiation is absorbed within the photosphere 
by molecules such as TiO, and so the stellar emis¬ 
sion is sub-blackbody at these wavelengths. As we 
will show later, the NUV/FUV plays a critical role 
in determining the steady state concentrations of 
GO and O 2 , rooted in the photochemistry outlined 
above. 


Returning to the model comparison, we note 
that the Tian et al. model also included an at¬ 
mospheric hydrogen budget. But when [Domagal-| 


Goldman et al. (2014) repeated these calculations. 


they got mixed results: these authors agreed with 
Hu et al. and disagreed with Tian et al. This out¬ 
come was somewhat surprising because the Hu et 
al. photochemical model was developed indepen¬ 
dently, but the Domagal-Goldman et al. and Tian 
et al. models were offshoots of the same Hast¬ 
ing group model. Thus, one might have expected 
the agreement/disagreement to be the other way 
around. Here, we compare our own (Kasting 
group) model with all three of these other mod¬ 
els in an attempt to identify the reasons for the 
discrepancies and to determine which predictions 
are robust and which are model-dependent. 

2. What Constitutes a False Positive for 
Life? 

Before proceeding further, we should step back 
and decide what exactly constitutes a false posi¬ 
tive for life. Oftentimes, this has been defined as 
detectable atmospheric O 2 or O 3 concentrations 
that are produced abiotically. But this definition 
is unsatisfying because the term ‘detectable’ de¬ 
pends critically on the technology used for detec¬ 
tion. Already, three different classes of telescopes 
have been advertised as having the ability to char¬ 
acterize Earth-like planets: i) JWST, ii) the TPF- 
type telescopes mentioned above, and iii) large 
ground-based telescopes. Each of these telescopes 
has different exoplanet characterization capabili¬ 
ties, and none of them may be able to measure 
atmospheric O 2 concentrations well below that of 
present Earth. But if any of these instruments 
were to succeed in identifying a promising, Earth¬ 
like planet, it is almost guaranteed that some later 
generation of astronomers would design a more ca¬ 
pable telescope to follow up on those observations. 
So, a more general way to think about this prob¬ 
lem is that any abiotic O 2 concentration that ex¬ 
ceeds measured, or inferred, biotically produced 
O 2 levels on Earth should be considered as a po¬ 
tential false positive. Whether it would be de¬ 
tectable by a first-generation TPF-type telescope 
is an interesting, but separate, question. 

We know very accurately what Earth’s O 2 con¬ 
centration is today - close to 21 percent by vol¬ 
ume (1 PAL). But we have good reason to believe 
that O 2 concentrations were much lower earlier 
in Earth’s history, including times during which 
oxygenic photosynthesis is believed to have been 
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operative. Oxygenic photosynthesis was invented 
by cyanobacteria during the Archean Eon, possi¬ 


bly as early as 3.0 billion years ago (Ga) (Crowe 
et al.|2013 1, although free O 2 did not begin to ac¬ 


cumulate in the atmosphere until around 2.5 Ga 
(the so-called Great Oxidation Event, or GOE). 
Some free O 2 was present during the Archean as 
a consequence of photolysis of CO 2 , but this O 2 
was largely confined to the stratosphere and was 
not appreciably more abundant than on the prebi- 
otic Earth (|Kasting et al.|1979| |Kasting & Catling 
20031). 


O 2 rose to appreciable concentrations in the 
lower atmosphere during the GOE, as evidenced 
by the disappearance of sulfur mass-independent 
fractionation in sedimentary rocks, along with 
other geologic O 2 indicators ( Holland|2006 1 . Dur¬ 
ing most of the ensuing Proterozoic Eon (2.5-0.54 
Ga), the atmosphere was oxygenated, but O 2 con¬ 
centrations are thought to have remained signif¬ 
icantly lower than today. Until recently, the hy¬ 
pothesized range of values for Proterozoic O 2 was 
0.01-0.5 PAL ( Kump|2008 1 . These values have al¬ 
ways been uncertain, however, and a recent study 
by Planavsky et al. (2014) suggests that mid- 


Proterozoic O 2 concentrations were at most 10“^ 
PAL. This estimate is uncertain, as it depends on 
the complex geochemical behavior of the element 
chromium, which has only recently been suggested 
as an O 2 indicator. However, this is currently our 
best estimate for Proterozoic O 2 , and we adopt it 
here as a lower bound for a biotically sustained O 2 
concentration. 

The Planavsky et al. results, if correct, im¬ 
ply that O 2 would have been difficult to de¬ 
tect directly during much of Earth’s history by 
a first-generation, optical, TPF-type telescope. 
Direct detection of the O 2 A band by such a 
device requires O 2 concentrations of ~0.01 PAL 
(Des Marais et al. 2002 Segura et al. 2003). In 


the visible wavelength region, then, the Protero¬ 
zoic Earth may be an example of an inhabited 
world with no detectable atmospheric biosigna¬ 
tures (Gockell 2014). However, a space telescope 


operating in the thermal-infrared could look for 
the ozone band at 9.6 /im, and that signal should 
be detectable for an O 2 level as low as 10“^ PAL 


(Leger et al. 1993 Segura et al. 2003). Alterna¬ 


be able to identify an analog to the Proterozoic 
Earth using a future space-based telescope oper¬ 
ating outside of the visible wavelength region. 

3. Atmospheric and Global Redox Bal¬ 
ance 

As mentioned earlier, many early attempts to 
calculate abiotic atmospheric O 2 concentrations 
did not properly account for redox balance. Redox 
balance is simply conservation of free electrons. 
Or, to say it another way, when one species is ox¬ 
idized (loses electrons), another species must be 
reduced (gain electrons). Both the atmosphere it¬ 
self and the combined atmosphere-ocean system 
must satisfy redox balance over time scales of tens 
of thousands to millions of years; otherwise, their 
compositions would change. The relevant time 
scale for redox balance in a weakly reduced early 
Earth-type atmosphere is ^30,000 yrs (Kasting 


2013). This is readily calculated by dividing the 


column depth of H 2 by the diffusion-limited es¬ 
cape rate of hydrogen. The relevant time scale 
for Earth’s 02 -rich modern atmosphere is about 


2 million years (Holland 2002). This is the reser¬ 


voir size of atmospheric O 2 divided by the rate 
of consumption of O 2 by surface weathering and 
oxidation of reduced volcanic gases. 

Atmospheric redox balance can be explicitly 
calculated in a photochemical model. Indeed, if 
the model is written self-consistently {i.e., con¬ 
serving mass and with stoichiometrically balanced 
chemical reactions), it must conserve redox bal¬ 
ance because free electrons can neither be gained 
nor lost. To track redox balance, we assign each 
species a redox coefficient based on the number 
of H 2 molecules that are generated or consumed 
when converting that species to a redox-neutral 
one. Following Kasting & Canfield (2012) and 
Kasting (2013), we define H 2 O, SO 2 , CO 2 , and 


tively, a large UVOIR telescope could potentially 
observe O 3 features in the UV. So, one might 


N 2 as neutral species for H- and 0-, S-, C-, and 
N-bearing gases, respectively. The redox coeffi¬ 
cient for each species can be found in Table in 
the Model Description (section]^. For example, 
H 2 S has a redox coefficient of -1-3, because H 2 S -I- 
2 H 2 O —>■ SO 2 -I- 3 H 2 , that is, one H 2 S molecule 
converted to SO 2 liberates 3 molecules of H 2 . An¬ 
other example would be H 2 O 2 , which has a redox 
coefficient of - 1 , as it consumes 1 molecule of H 2 
when converted to water vapor (H 2 O 2 -I- H 2 —>■ 


6 



















































2 H 2 O). Under these assumptions, the atmospheric 
redox balance can be written as 

^ out (Red) + 

^rain (Ox) = ^esc(H 2 ) + ^rain (Red) 

( 1 ) 

Here, ^out (Red) is the outgassing rate of reduced 
species in all chemical forms, weighted by their 
stoichiometric coefficient for H 2 production 

^out (Red) = ^out(H2) + ^ out (CO) 

+i<^out(CH4) + 3<^out(H2S) + ... ( 2 ) 

^esc(H 2 ) is the hydrogen escape rate to space, as¬ 
sumed to be given by the diffusion-limited flux 


Unlike atmospheric redox balance, a photo¬ 
chemical model will not compute global redox bal¬ 
ance automatically because it does not include an 
ocean. But global redox balance can be imposed 
by making certain assumptions about the magni¬ 
tude of the terms in eq. |^and then adjusting the 
lower boundary conditions of the photochemical 
model so that the retained terms are in balance. 
On an abiotic planet, oxidative weathering and all 
sediment burial terms are expected to be small. If 
we neglect these terms, then eq. [^simplifies to 

^out(Red) = ^esc(H2) (5) 


^esM) « 2.5 X 10^^fT(H2) 


(3) 


and fT(H 2 ) (= 0.5f(H) + f(H 2 ) + /(i?20) + 
2 f(CH 4 ) -|- ...) is the total volume mixing ratio 
of hydrogen in all its chemical forms. Meanwhile, 
^rain(Ox) and ^rain(Red) are the combined rain- 
out and surface deposition rates of oxidants and 
reductants, respectively. We treat rainout and 
surface deposition together because both processes 
are assumed to result in transfer of gases from the 
atmosphere to the ocean. 

Global redox balance refers to redox balance 
within the combined atmosphere-ocean system. In 


the notation of Kasting (2013), this can be ex¬ 
pressed as 

*^out(Red) + $OW + *^burial(CaS04) + 

^burial (Re^O4) = *^esc(dd 2 ) “t" burial (CH2O) 

+ 5<^bur^al(FeS2) (4) 

H 2 sources are on the left-hand side of this equa¬ 
tion; H 2 sinks are on the right. Here, ^ow 
represents oxidative weathering of the continents 
and seafloor and ^buriai(i) is the burial rate of 
species i in sediments. Important species are gyp¬ 
sum (CaS 04 ), magnetite (Fe 304 ), organic mat¬ 
ter (CH 2 O), and pyrite (FeS 2 ). We have written 
them in this manner to show the stoichiometry 


explicitly. Other authors (e.g.. Catling & Claire 


20051 have used different notation. Many authors 


e.g., Holland 2002 2009) compute the global re¬ 


dox budget in units of O 2 , rather than H 2 , and 
they adopt sulfate (SO 4 ) as the reference oxida¬ 
tion state for sulfur (although frequently these as¬ 
sumptions are not stated explicitly). These latter 
parameters are good choices for the 02 -rich mod¬ 
ern atmosphere but are less well suited for I 0 W-O 2 
atmospheres. 


i.e., volcanic outgassing of reduced gases must be 
balanced by escape to space. By combining this 
expression with eq. we see that our model must 
also satisfy the equation 


(f) 


i(Ox) — ^j-oin(Red) 


( 6 ) 


that is, rainout plus surface deposition of oxidants 
must equal rainout plus surface deposition of re¬ 
ductants. Another way to think of this is that, if 
consumption of oxidants by the crust and burial 
of reductants in sediments are neglected, oceanic 
redox balance requires that the inputs of oxidants 
and reductants from the atmosphere must be in 
balance. We will return to this thought below 
because it is key to understanding the modeling 
approach taken here, as well as those adopted by 


Tian et al. 

(2014 

) and 

Domagal-Goldman et al. 


(2014). Until those articles appeared, most pa¬ 


pers regarding false positives (including some writ¬ 
ten by us, e.g., Segura et al. (2007)), have simply 


overlooked the concept of global redox balance. 

4. Model Description 

4.1. Model parameters and atmospheric 
chemistry 


To compare with other recent false positive cal¬ 
culations, we used a 1-D, horizontally averaged, 
photochemical model that shares its heritage with 
that of Segura et al. (2007), Tian et al. (2014), 


and Domagal-Goldman et al. (2014). The model 
has 30 long-lived species (Table [^, 15 short-lived 
species, and one aerosol species (sulfate). N 2 is 
treated as inert and is given a constant vertical 
mixing ratio set by the abundance of other gases 
so as to ensure a 1-bar surface pressure. This 
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keeps us out of the parameter space in which the 
background pressure is low and water loss is rapid 
(Wordsworth & Pierrehumbert 2014). Chemical 


steady state is calculated via 215 reactions, includ¬ 
ing the photolysis of key atmospheric constituents. 
The full reaction list is given in the Appendix. 

We did not attempt to calculate self-consistent 
temperature profiles for different model atmo¬ 
spheres. Rather, we assumed temperatures com¬ 
parable to those found in the models with which 
we wish to compare (e.g., |Hu et al.||2012l |Tian 


et al.|[20l4 ). Similarly, we employed two different 
eddy diffusion profiles: the first from Hu et al. 
for our 90% CO 2 cases, and the second from Tian 
et al. for our 5% CO 2 cases. The eddy diffusion 


profile used by Tian et al. (2014) results in a low¬ 


ered homopause (located around 75 km), which 
accounts for the decreasing CO 2 concentrations in 
the upper atmosphere of their model. 

4.2. Stellar spectra 

Some of our calculations involve Earth-like ex¬ 
oplanets orbiting F, K, and M stars. In such 
cases, we have adopted the methodology of |Se-| 


gura et al. 

( 

2005 

1 to manipulate the stellar spec- 

tra. From 

Segura et al. 

(2005), we normalize and 


scale the spectrum of AD Leonis (AD Leo), a well- 
studied, active M3.5Ve star (a star that periodi¬ 
cally/episodically flares from a quiescent state), 
as well as e Eridani, a KIV star, and a Bootis, 
an F2V star. Additionally, we have prepared the 
spectra of the M star GJ 876 from the short¬ 
wave observations of the MUSCLES observations 
( France et al.||2013 ) coupled to a NextGen model 
for a 3,200 K model star. These spectra, along 
with the absorption cross-sections for important 
atmospheric constituents, are shown in Fig. [l] 
Additional details about each star can be found 
in Table 121 

4.3. Spectral models 

To generate synthetic direct imaging spectra 
of planets with the atmospheres modeled in this 
study, we used the Spectral Mapping Atmospheric 
Radiative Transfer Model (SMART), developed by 
D. Crisp (|Meadows & Crisp||1996 [Crisp 1997). 


SMART has been well-validated against spacecraft 


observations of the Earth (Robinson et al. 2011 


2014) and Venus (Arney et al. 2014). An enhanced 


version of SMART that includes the effects of re¬ 
fraction is used to generate synthetic transit trans¬ 
mission observations of our model atmospheres 
(Misra et al. 2014a|b ). SMART’S transit trans¬ 
mission capability has been validated against both 
solar occultation soundings of the Earth’s atmo¬ 
sphere ( Gunson et al.||1990 l and Earthshine spec¬ 
tral measurements ( Palle et al.|2009 ). The absorp¬ 
tion cross-sections in the visible to near-infrared 
wavelength range (A >0.4 /xm) for each spectrally 
active gas are calculated by the SMART routine 
Line-By-Line ABsorption Coefficients (LBLABC) 
using the HITRAN 2012 line list database ( [Roth- 


|man et 2013). For UV wavelengths (A <0.4 
/xm), we use the absorption cross sections in Fig. 

To calculate the synthetic transit transmis¬ 
sion spectra, we use an index of refraction of n 
= 1.00031 at STP (consistent with an atmospheric 
mixture of ^95% N 2 and ~5% CO 2 ) and a planet- 
star impact parameter of 6 = 0 (where the planet 
is assumed to transit directly through the center 
of the star). For the synthetic direct imaging cal¬ 
culations, we assume the surface is a Lambertian 


ocean (McLinden et al. 1997) and the planet is 
observed with a solar zenith angle of 60°, approxi¬ 
mating a disk-averaged observation at quadrature 
phase ( Segura et al.|2005 ). We used the stellar pa- 
rameters given in Table The spectral resolution 
for each case is A=1 cm^, which is approximately 
AA= 5.8x10“® /xm at A=0.76 /xm. 

4.4. Deposition velocities and ocean chem¬ 
istry 

Soluble long-lived species are rained out of the 
troposphere using the parameterization of |Giorgi| 
& Chameides (1985). In some cases, when the 


ocean is assumed to be saturated with a particu¬ 
lar gas, rainout is turned off because it would nec¬ 
essarily be balanced by an upward flux from the 
ocean surface. Long-lived species may also have 
a non-zero deposition flux if the ocean is not sat¬ 
urated with that gas. Some detail is given here, 
as this process can directly affect the atmospheric 
and global redox balance. Indeed, as we will show, 
for planets orbiting M stars, the atmospheric O 2 
abundance depends critically on the deposition ve¬ 
locities for O 2 and CO, and these quantities, in 
turn, depend on the chemistry of these gases in 
solution. 

For minor, highly reactive gases and radicals. 







































































Table 1 

Chemical species 


Chemical 

Deposition 

Redox 

formula 

velocities (cm s“^) 

coefficient 

0 

1 . 

-1 

O 2 

(0-1.4x10-4) 

-2 

H 2 O 

fixed^ 

0 

H 

1 . 

0.5 

OH 

1 . 

-0.5 

HO 2 

1 . 

-1.5 

H 2 O 2 

0.5 

-1 

H 2 

variable^ 

1 

CO 2 

fixed‘s 

0 

CO 

(0-1.2x10-4) 

1 

HCO 

1 . 

1.5 

H 2 CO 

0.1 

2 

CH 4 

0 . 

4 

CH 3 

1 . 

3.5 

C 2 H 6 

1 . 0 x 10 -'^ 

7 

NO 

3.0x10-4 

-1 

NO 2 

3.0x10-4 

-2 

HNO 

1 . 

-0.5 

H 2 S 

0.015 

3 

HS 

3.0x10-4 

2.5 

S 

1 . 

2 

HSO 

1 . 

1.5 

SO 

3.0x10-4 

1 

SO 2 

1 . 

0 

NH 3 

0 . 

1.5 

NH 2 

0 . 

1 

N 

0 . 

0 

N 2 H 4 

0 . 

2 

N 2 H 3 

0 . 

1.5 

H 2 SO 4 

0.2 

-1 


Note. —^ H 2 O volume mixing ratio is fixed; ^ 
H 2 is used to balance the global redox budget - 
see the text for a description; The CO 2 mixing 
ratio is fixed between 0.05 and 0.9. 
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Table 2 

Stellar properties 


Star 

Spectral 

type 

Te// 

(K) 

Luminosity 

L0 

Stellar 
Radius [R©] 

Distance 

(pc) 

Semi-major 
axis (AU) 

Ly-o at 1 AU 
(ergs cm~^ 

a Bootis 

F4Vi 

6435 I 

3 .I 5 I 

1.431^ 

15.8® 

2.31 

79 . 7 ® 

Sun 

G2V 

5780 

1 

1 

- 

1.3 

8.17 

e Eridani 

KIV'* 

5039® 

0.32® 

0.735® 

3.2“ 

0.76 

47.1"^ 

AD Leonis 

M3.5Ve4 

3390® 

0.023® 

0.39^ 

4.7® 

0.21 

7.05® 

GJ 876 

M4V® 

3129® 

0.012® 

0.38® 

00 

0.15 

0.414® 


Note.— ^ iBoyaiian et al.| (|2013||; ^ iJanson et al.|(|2013 

; ^iLandsman & Simon] (ll993|l: ^]Linskv et al.] | 

Baines &: Armstrong]||2012|1; ° jRojas- Ayala et al.|(|2012||; ‘ 

Reiners et al.|l|2009|l; ®|von Braun et al.|l|2014||; 

et al.|(|2013f. 


France| 


we assume constant deposition velocities ranging 
from 0.01-1 cm s“^ derived from the modern Earth 
literature (e.g., Liss & Slater 19741. These de¬ 


surface ocean, and C (= 6.02 x 10^° L mol ^ cm 


position velocities are limited by the rate of gas 
transfer within the atmospheric boundary layer. 
For longer-lived constituents such as O 2 and CO, 
we use a different approach. Because we are in¬ 
terested in gas transfer to the ocean, rather than 
surface weathering, we assume that the deposition 
velocity of such species is controlled by its ‘piston 
velocity’ within the oceanic boundary layer, 


[cm s ^]) (|Broecker & Peng|l982 ^ 

Charecha et al. 

2005 Domagal-Goldman et al.|201^ 

Technically, 


the oceanic boundary layer resistance adds to the 
atmospheric resistance, but because the oceanic 
resistance is much larger, we can simply ignore the 
atmospheric boundary layer for long-lived species. 
The piston velocity, in turn, is determined by the 
vertical gradient in species concentration within 
a thin film at the top of the ocean and by its 
molecular diffusivity, ^diff (cm^ s“^), in solution. 
Mathematically, Vp = Kdiff/zfum, where ZfUm is 
the (empirically determined) thickness of the film, 
~40 fim. 

To give a specific example, the depositional flux 
of carbon monoxide into the ocean is calculated 
from the expression: 

^{CO) = Vdep{CO)nco = Vp {acoPco — [CO]s)-C, 

(7) 

Here, aco (= 1-0x10 ^ mol L ^ at 25°C) is the 
Henry’s Law coefficient for CO, pco is the partial 
pressure of CO at the bottom of the atmosphere, 
[CO]s is the assumed concentration of CO in the 


The piston velocity, 
-1 


is a unit conversion factor. 

Vp, for CO is 4.8x10“^ cm s“^ (Kharecha et al. 
|2005[ ) . The deposition velocity (the values seen in 
Table[^ is then given by VdepiCO) = ^(CO)/nco, 
where nco (=2.5x10^® cm“^ x fco) is the num¬ 
ber density of CO at the lowermost grid point. 
Each calculated deposition velocity bears with it 
an assumption about the dissolved concentration 
of that species in the ocean. An upper limit on 
species deposition velocity is obtained if the dis¬ 
solved concentration is equal to zero. For CO, at 
25°C, this value is ~1.2xl0“'^ cm s“^. O 2 has 
a similar solubility and diffusivity, so under the 
same conditions its maximum deposition velocity 
is similar (~1.4xl0“^ cm s“^). 

For species whose concentrations in the ocean 
are non-zero, calculating the deposition velocity is 
more complicated. Kharecha et al.l (20051 calcu¬ 


lated abiotic deposition velocities for CO using a 
2-box ocean model (surface and deep), along with 
a limited set of aqueous reactions. They estimated 
a value of 10“®-10“® cm s“^, assuming that the 
largest sink for CO in the ocean was the hydration 
of CO to formate (CO -I- OH“ —HCOO“). Here, 
we expand the aqueous chemistry of Kharecha et 


al. to include reactions between dissolved O 2 , 
and formate: 

CO, 

I/ 2 O 2 + CO 

—t 

CO2 

(8) 

CO + H2O 

—t 

CO2 H2 

(9) 

I/2O2 + HCOO- 

—t 

OH- + CO2 

(10) 


These reactions could potentially take place in the 
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ocean or as CO/02-bearing water flows through 
midocean ridges. The first two reactions have been 


studied in research on catalysis (see Hibbitts & 

M 


:lesia 


(2015) for a review); they represent the di¬ 
rect oxidation of CO (eq. also included in Fig. 

and the water-gas shift reaction (eq. [^. An 
extensive literature survey has not revealed any 
studies that test these reactions in the absence of 
catalysts, and we will come back to this issue in 
the Discussion section. The third suggested reac¬ 
tion (eq. 101 is the direct oxidation of formate, 


a process that can be inferred from radiochemical 
studies (e.g., Balkas et al.|1966 | and investigations 
of supercritical water reactions (Helling & Tester 


1988). A moderate amount of atmosphere attenu¬ 


ates X-rays enough to preclude the radiochemical 


pathway, even for more active stars (Cnossen et al. 


2007), but this reaction may still occur slowly in 


the absence of X-rays. These reactions have been 
incorporated into the model of Kharecha et ah, as 
outlined in the Appendix, and will be used to con¬ 
strain the deposition velocities for CO and O 2 in 
the Results section. 

4.5. Balancing the global redox budget 

As discussed in Section a photochemical 
model will not satisfy global redox balance auto¬ 
matically; however, one can impose global redox 
balance by imposing appropriate boundary condi¬ 
tions. In this study, we are primarily interested 
in calculating upper limits on abiotic O 2 . Hence, 
in eq. we neglect the O 2 sinks from oxida¬ 
tive weathering ($oiv) and from burial of gypsum 
(‘hbMriai(CaS 04 )). In many of our calculations we 
neglect the O 2 sink due to ferrous iron oxidation 
(d>b«riai(Fe 304 )), as well, even though this term 
would arguably be important on an Earth-like ex¬ 
oplanet. We also neglect burial of organic car¬ 
bon ($b„„a/(CH 20 )) and pyrite (4>6„r*a/(FeS2)). 
These O 2 sources are important on the modern 
Earth, but they are both primarily biological. Or¬ 
ganic carbon is formed directly by organisms, and 
pyrite is formed by biologically catalyzed reaction 
of this organic carbon with dissolved sulfate in the 
ocean. Oceanic sulfate levels should be low on an 
abiotic planet, as they were on the Archean Earth 


prior to the rise of atmospheric O 2 (Canfield et al. 


2000 ). 


Under these assumptions, the global redox bud¬ 
get can be represented by eq. Rainout of ox¬ 


idants from the atmosphere must equal the rain¬ 
out of reductants. To balance this equation, we 
have implemented a scheme similar to that used 
by |Tian et'n(] ( |2014| ). If ^rain{Ox) < ^rain{Red), 
which is the case for the overwhelming majority 
of our simulations, the excess hydrogen flowing 
into the ocean is returned to the atmosphere as 
H 2 . We accomplish this by imposing an upward 
flux of H 2 at the lower boundary. H 2 is also out- 
gassed from volcanoes, but we simulate this pro¬ 
cess (and other volcanic outgassing rates) by in¬ 
troducing fictitious chemical source terms within 
the model troposphere. This procedure allows us 
to use the lower boundary condition on H 2 to bal¬ 
ance the global redox budget. In the case where 
^rain(Ox) > ^rain{R^d), we impose an increased 
depositional velocity for H 2 . This latter case is 
rarely observed, though, because the weakly re¬ 
duced atmospheres studied here consistently pro¬ 
duce more soluble reductants than oxidants. All 
of the results presented here have balanced global 
redox budgets. 

In some cases, we have assumed a non-zero fer¬ 
rous iron oxidation term, (Fe 304 ), in eq. 

(the reaction of ferrous iron with O 2 is also in¬ 
cluded in Fig. [^. Oxidation of dissolved fer¬ 
rous iron is a quick and readily available O 2 sink, 
as has been recognized for many years in stud¬ 


ies of the early Earth (e.g., Cairns-Smith 1978 


Braterman et al. 1983 Kasting et al. 1984). Bi¬ 


ology is not needed to accomplish this process, as 
it can be catalyzed by near-UV sunlight (ibid.). 
The deposition of banded iron-formations (BIFs) 
provides evidence that ferrous iron was abundant 
in the Archean oceans (Holland 1973), and this 
would probably be true on abiotic exoplanets, as 
well. That said, most BIF deposition on Earth 
was probably triggered biologically. 


5. Results 

5.1. Planets around G stars 

5.1.1. Planets with high H 2 outgassing rates 

We begin by pointing out that our model pro¬ 
duces very little ground-level O 2 for typical abiotic 
Earth conditions, regardless of the CO 2 concentra¬ 
tions in the atmosphere. Here, a typical abiotic 
Earth would resemble the prebiotic Earth, which 
is thought to have had more CO 2 than the mod- 
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ern Earth. The adjustment would be caused by 
the negative feedback between CO 2 and climate 
(Walker et al. 1981). The prebiotic Earth might 
have had ~80% N 2 and enough CO 2 (^ 20 %) to 
provide clement conditions. Assumed volcanic 
outgassing rates for the early Earth are ~3xl0® 
cm“^ s“^ for SO 2 , ~3xl0® cm“^ s“^ for H 2 S, 
~3xl0® cm“^ s“^ for CH 4 , and ~10^° cm“^ s“^ 
for H 2 (Segura et al. 2007). Under these con¬ 
ditions our model surface O 2 concentrations are 
^1x10“^^ PAL, similar to those shown in Fig. 4a 
of Segura et al. (2007). Outgassing rates within 
an order of magnitude higher or lower than those 
quoted here seem reasonable for an Earth-like 
planet throughout most of its lifetime, and pro¬ 
duce similar (very low) surface O 2 concentrations. 
Here, an ‘Earth-like’ planet is defined as one that 
has a reduced interior, active volcanism, and a 
volatile inventory similar to the Earth’s. 


5.1.2. Planets with low outgassing rates 


We next changed the parameters in the model 
to try to duplicate false positives that have been 
published in the recent literature. The first case 
we looked at was a high-C 02 atmosphere (90% 
CO 2 , 10% N 2 ) on a planet orbiting a G star like 
our Sun. Hu et al. (2012) simulated such a planet 
with three different H 2 outgassing rates: 3x10^°, 
3x10®, and 0 cm“® s“^. Their zero H 2 outgassing 
case still included a small outgassing flux of H 2 S, 
as discussed further below. The assumed depo¬ 
sition velocities for CO and O 2 were 1x10“® cm 
s“^ and zero, respectively. Their results are shown 
in the two left-hand panels of Fig. Surface 
O 2 was low in the high-outgassing case, but it 
increased markedly as volcanic outgassing rates 
decreased. Ultimately, their zero-outgassing case 
yielded a well-mixed O 2 profile with a mixing ra¬ 
tio of 1.3x 10“®. As pointed out earlier, this is just 
slightly above the inferred O 2 mixing ratio for the 
Proterozoic Earth, so this constitutes a false pos¬ 
itive for life by our definition. 


Using the exact same boundary conditions on 
CO and O 2 , we do not generate high surface O 2 
concentrations, for reasons that remain unclear 
(see discussion below). These boundary conditions 
do not yield a balanced global redox budget, and 
so we have chosen not to show them. However, 
these results differed only slightly from the same 
calculation with a fully balanced global redox bud¬ 


get and we have shown these latter calculations in 
the right-hand panels of Fig. The global redox 
budget for these three balanced cases is shown in 
Table The surface O 2 mixing ratio in our model 
never rises above 10 “^® in any of our simulations; 
instead, ground-level O 2 mixing ratios are always 
low, as they are in the high-outgassing case. 

Why would Hu et al. predict 10 orders of mag¬ 
nitude more ground-level O 2 than we do for the 
low-out gassing planet? The answer is not clear. 
These two photochemical models were developed 
independently, and both contain many different 
parameters. The Hu et al. model has a more 
extensive reaction set and has different rate con¬ 
stants for some reactions, especially those involved 
in production of CH 4 . It is not obvious, though, 
why this should cause any significant difference 
in the results. A model intercomparison study is 
planned for the coming year to try to resolve this 
discrepancy. 

We can make two observations, though, that go 
beyond these technical differences in the computer 
models. First, low-outgassing. Earth-sized plan¬ 
ets are probably not physically realistic. Mars is 
a low-outgassing planet today, but that is because 
it is only one-tenth of Earth’s mass and is vol¬ 
canically inactive. A true Earth analog would be 
unlikely to be this dead unless it was extremely 
old. Second, a planet that was totally volcanically 
inactive, so that it released zero reduced gases at 
its surface, could have a much larger O 2 concen¬ 
tration than the one found by Hu et al. According 
to eq. if we neglect surface oxidation and sedi¬ 
ment burial terms for such a planet, the volcanic 
outgassing rate must equal the hydrogen loss rate 
to space. A zero-outgassing planet must therefore 
have a zero hydrogen escape rate if its atmosphere 
is in redox steady state. But Earth-like planets, 
by definition, are considered to have liquid water 
on their surfaces and some water vapor in their at¬ 
mospheres. This alone is enough to create a finite 
hydrogen escape rate, leaving behind the oxygen 
from the water. Thus, O 2 should actually accu¬ 
mulate indefinitely on such a planet. One reason 
it does not in the Hu et al. model, and in our own 
calculation shown in Fig. is that both calcula¬ 
tions assumed a finite outgassing rate of H 2 S, even 
while other outgassing rates were set to zero. H 2 S 
can be converted photochemically to H 2 according 
to the stoichiometry: H 2 S -I- 2 H 2 O —SO 2 -I- 3 
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Baseline 

Reduced 

H2S only 

H2 Escape 

3.21x101® 

5.10x10® 

9.00x10® 

Red. Outgassing 

3.21x101® 

5.10x10® 

9.00x10® 

^ OX 

3.57x10® 

6.98x10® 

1.17x10® 

SO4 (aerosol) 

3.57x10® 

6.90x10® 

1 .11x10® 

H2SO4 

0 

7x10® 

2 x101' 

H2O2 

0 

9x10® 

3x101' 

^red 

4.58x10® 

3.30x10® 

1.67x10® 

HSO 

3.53x10® 

2 .68x10® 

1.26x10® 

H2S 

5.1x10® 

3.9x10® 

2 .6x10® 

H2CO 

3.5x10® 

1.3x10® 

8 x101' 

HCO 

1.3x10® 

9x101' 

4x101' 

CO 

6 x101' 

IxlOi' 

3x10® 

Ee+“'' Elux 

0 

0 

0 

H2 Return 

4.22x10® 

2.60x10® 

4.96x10® 

H2 Balance 

IxlQi^ 

3x10^ 

1 x10^ 

O2 column depth 

2 .1x101® 

3.3x101® 

4.5x101® 

O3 column depth 

3.9x101'! 

1.5x101® 

4.9x101® 


H 2 . So, neither model actually represents a case 
of zero outgassing of total hydrogen. 

The atmospheric redox budget in both models 
is influenced by rainout and surface deposition of 
oxidants and reductants at the atmosphere-ocean 
interface. These processes are particularly impor¬ 
tant in the Hu et al. model because of the large 
amount of O2 present in the lower atmosphere. 
To first order, surface deposition of the oxidants 
O3, HO2, and H2O2 in their model is balanced 
by deposition of the reductant CO. (They com¬ 
pute this balance in units of H, instead of H2, but 
the implications are the same.) So, if one were to 
rethink their model in terms of global redox bal¬ 
ance, this implies that CO must be reacting with 
these oxidants in solution within the ocean. That, 
by itself, is not physically implausible. But O2 
in their model is given a zero deposition velocity 
at the surface, for reasons that are not explained. 
What if O2 itself reacts with CO (or formate) in 
solution? Its deposition velocity would then be 
non-zero, and its ground-level mixing ratio would 
presumably be much lower. As we discuss further 
below, the rate of aqueous reaction of CO and O2 
is a key parameter in deciding whether false posi¬ 
tives are physically realistic. Our own model does 
not generate a false positive in this situation, so we 
will not spend further time analyzing it. For now, 
we simply use this example to illustrate why pro¬ 
cesses happening below the bottom of one’s pho¬ 
tochemical model must be taken into account. 


Table 3: Global redox budget for the three cases in 
the right panels of Fig. Baseline: 3x10^° cm“^ 
s“^ H2, 3x10® cm“^ s“^ CH4, 3x10® cm“^ s“^ 
H2S; reduced: 3x10® cm“^ s“^ H2, 3x10® cm“^ 
s“^ CH4, 3x10® cm“® s“^ H2S; H2S only: 3x10® 
cm“^ s“^ H 2 S. All fluxes are in units of cm“® s“^. 


5.2. Planets around M stars 

A second possible false positive could occur on 
planets orbiting M stars. We already dismissed 
such planets as being habitable in the Introduc¬ 
tion, because three different studies have predicted 
that these planets would lose most or all of their 
water during the pre-main-sequence phase. That 
type of false positive, we argued, would be easy 
to identify because of the lack of H2O absorption. 
But let us imagine an M-star planet that either 
started off with much more water than Earth (e.g., 
and sufficient sinks to remove 
up from its pre-main sequence 
time, or a planet that migrated inward to the hab¬ 
itable zone at some time after the parent star had 
evolved onto the main sequence. Could such a 
planet accumulate a signihcant amount of O2? 

This type of planet has been modeled by |Tian| 


huger et al. 2015 


the O2 that built 
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et al. (2014) and Domagal-Goldman et al. (2014), 


with conflicting results. Both groups studied a 
planet orbiting the M star, GJ 876. Tian et al. 
predict that a planet with 5% GO 2 in its atmo¬ 
sphere and Earth-like volcanic outgassing rates 
could reach an O 2 surface mixing ratio of a little 
more than 2x10“^ (see their Fig. 3b). This, by 
our definition, is a false positive. A similar planet 
orbiting the Sun would have a surface O 2 mixing 
ratio of more like 10“^^. (This value is not actu¬ 
ally shown on their Fig. 3a, but it can be inferred 
from our own model, which gives similar results.) 
The Tian et al. results for the M-star planet are 
strongly influenced by their lower boundary condi¬ 
tions, as discussed further below. Guriously, when 
Domagal-Goldman et al. model the same M-star 
planet, using identical boundary conditions, they 
find 10^ times less O 2 . When they use their own 
lower boundary conditions to find the maximum 
abiotic O 2 production, O 2 increases by less than a 
factor of 10, leaving it still 10^ times lower than the 
Tian et al. model. Glearly, lower boundary condi¬ 
tions are important, but other detailed differences 
between these two different versions of the same 
parent photochemical model may also be influenc¬ 
ing the results. 


We have qualitatively duplicated the results of 
Tian et al. in our Fig. Our model predicts 
surface CO and O 2 mixing ratios of ^7x10“^ and 
~3xl0“^. The slight difference in CO and O 2 
mixing ratios compared to their results is a func¬ 
tion of the assumed distance to the parent star: 
because the abiotic Earth modeled by Tian et al. 
is further from GJ 876, photolysis rates are lower, 
resulting in less CO and O 2 . We note also that 
Tian et al. have either an implied methane flux or 
a different set of rates than those reported here; 
the calculated methane mixing ratio in our model 
(with no methane outgassing) is about 2 orders of 
magnitude smaller than theirs for the solar case 
and is completely negligible for the GJ 876 case 
(off the scale to the left for both stars in Fig. |^. 
Aside from this, our calculated CO and O 2 surface 
mixing ratios are quite close to theirs for both the 
G-star planet and the M-star planet. This agree¬ 
ment may be caused, in part, by the fact that 
we balanced the global redox budget in the same 
way that they did - by adding an upward flux of 
H 2 from the ocean when rainout of reductants ex¬ 
ceeded rainout of oxidants. 


That said, the agreement between our calcula¬ 
tions and those of Tian et al. does not imply that 
the results are realistic. The calculated O 2 mixing 
ratios for the M-star planet are strongly influenced 
by the assumed lower boundary conditions which, 
in this case, are somewhat arbitrary. Following 
Tian et ah, we used an Earth-like rainout rate 
and deposition velocities of 10“® cm s“^ for both 
CO and O 2 . The value for CO is approximately 
the logarithmic mean of its maximum deposition 
velocity of 1.2x10“^ cm s“^ and the abiotic de¬ 
position velocity of ~10“® cm s“^ calculated by 


Kharecha et al. 

(2005 

) (see Sect. 

4.1 


not a particularly strong justification for choos¬ 
ing this value. The O 2 deposition velocity was 
set to the same value, again for no clearly iden¬ 
tifiable reason. These deposition velocity choices 
have little effect on the G-star calculation because 
the surface concentrations of both O 2 and CO 
are kept low by photochemistry (Fig. |^. But 
for the M-star case, they give rise to large down¬ 
ward fluxes of both species. In our model, the sur¬ 
face O 2 mixing ratio is 3.2x10“^ and the surface 
CO mixing ratio is 6.6x10“^ (Fig. |^. The to¬ 
tal number density in the lowermost model layer 
is ~2.5xl0^®, so the downward deposition flux, 
= V^depf^ntot, IS 8.0x10^° Cm-^ S"! 
for O 2 and 1.7x10^^ cm“^ s“^ for CO. The flux 
of CO is almost exactly equal to twice that of O 2 , 
which is not surprising, given that they are pro¬ 
duced in this same ratio by CO 2 photolysis fol¬ 
lowed by reaction of O with itself to produce O 2 . 
This also implies that these gases must be reacting 
within the ocean according to: 2 CO -I- O 2 —2 
CO 2 . As can be seen from Table|^ the CO and O 2 
deposition fluxes dominate the atmospheric redox 
budget in this case, so the answer that one gets is 
very much a function of the assumed lower bound¬ 
ary conditions on O 2 and CO. 

Let us consider what these chosen deposition 
velocities imply about ocean chemistry. Accord¬ 
ing to eq. the fact that the chosen deposition 
velocity is ~1 percent of the maximum deposition 
velocity means that the ocean must be 99 percent 
saturated with both CO and O 2 . That, in turn, 
implies that the rate of reaction of these species is 
very slow. If the rate of reaction of CO and O 2 in 
solution was fast, then their deposition velocities 
should be ^10“^ cm s“^, and their surface mixing 
ratios would drop by more than a factor of 200. 
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The reaction rate between dissolved CO and O 2 is 
discussed further in Section r5.2.21 


5.2.1. ‘Worst case’ scenarios leading to signifi¬ 
cant false positives 

One complication in the methodology being em¬ 
ployed is that CO and O 2 are transferred from 
the atmosphere to the ocean by two different pro¬ 
cesses - rainout and surface deposition. For the 
G-star planet, neither species is produced in ap¬ 
preciable amounts, so both of these fluxes are neg¬ 
ligible. For the M-star planet, both gases are much 
more abundant, and one needs to think carefully 
about what this assumption implies. For exam¬ 
ple, a worst-case scenario {i.e., a likely false posi¬ 
tive for life) would be if the ocean was completely 
saturated in O 2 . Following eq. this would im¬ 
ply a zero deposition velocity for O 2 . However, 
the net rainout rate should be zero in this case as 
well, because any O 2 entering the ocean by way 
of raindrops would necessarily be balanced by an 
upward flux of O 2 from the ocean surface. So, 
the real upper limit on abiotic O 2 comes from set¬ 
ting both rainout and surface deposition equal to 
zero for O 2 . We have illustrated such a case in 
Fig. [^for the Sun, GJ 876 (an M star), as well 
as for a Bootis (an F star), e Eridani (a K star), 
and AD Leonis (another M star). Table [^details 
the global hydrogen budget for each case. All of 
these cases (save for the planet around AD Leo, 
which is similar to the one around GJ 876) have 
been used to generate spectra (Fig. [^. For these 
cases, we kept the deposition velocity for GO at 
10 “® cm s“^, which is equivalent to assuming that 
formate reacts directly with oxidants in the ocean 
but that GO does not; thus, the rate-limiting step 
in GO removal from the surface ocean is hydration 
to formate. 



Sun 

GJ 876 

H 2 Escape 

l.dxlQi'J 

1.4x101“ 

H 2 Outgassing 


nr™ 

^OX 

1.52x10*^ 

2.36x1011 

O 2 

0 

2.34x1011 

HO 2 

2 x 10 ® 

9x10“ 

H 2 O 2 

1.50x10® 

2 x 10 “ 

^red 

6.45x10“ 

2.38x1011 

GO 

4.0x10® 

2.38x1011 

HGO 

2.3x10® 

0 

H 2 GO 

5.82x10“ 

0 

Ee++ Elux 

4x10“ 

4x10“ 

H 2 Return 

1.03x101“ 

6.16x10“ 

H 2 Balance 

3x10“ 

2 x 10 ^ 

O 2 column depth 

1 . 6 x 101 “ 

9.3x10^^ 

O 3 column depth 

5.4x1014 

8.4x1011’ 


Table 4: Global redox budget for Fig. All fluxes 
are in units of cm“^ s“^. 
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a Bootis 

Sun 

e Eridani 

AD Leo 

GJ 876 

GJ 876“ 

H2 Escape 


lo™ 

To™ 

To™ 

lo™ 

lo™ 

H2 Outgassing 

10 ™ 

rF 

To™ 

IcP 

10 ™ 

10 ™ 

^ OX 

7.01x10“ 

2.52x10“ 

7.30x10“ 

9.88x10“ 

1.43x10“ 

4.76x10““ 

O2 

0 

0 

0 . 

0 

0 

2 .12x10““ 

HO2 

3x10“ 

2 x10“ 

2 .1x10® 

1 x10'^ 

5x10“ 

3x10“' 

H2O2 

3.37x10“ 

2.50x10® 

7.09x10“ 

9.87x10“ 

1.42x10“ 

1 x10“ 

^red 

6.93x10™ 

5.72x10“ 

1.06x10““ 

2 .01x10““ 

2.92x10“ 

4.77x10““ 

CO 

2 x10“ 

3x10“ 

1.06x10““ 

2 .01x10““ 

2.92x10“ 

4.77x10““ 

HCO 

1 .2x10“ 

2 .2x10® 

0 

0 

0 

0 

H2CO 

6.81x10™ 

5.49x10“ 

0 

0 

0 

0 

H2 Return 

6.93x10™ 

5.47x10“ 

3.30x10“ 

1.03x10““ 

1.49x10“ 

1 .00x10“ 

H2 Balance 

-5x10^ 

5x10“ 

4x10“ 

-2x10® 

-1x10® 

8 x10“ 

O2 column depth 

6.4x10™ 

1 .8x10™ 

1 .0x10““ 

2.5x10“® 

1.4x10““ 

1.5x10““ 

O3 column depth 

8.4x10™ 

6 .8x10™ 

4.3x10““ 

2.9x10“® 

1.3x10“® 

1 .1x10““ 


Table 5: Global redox budget for Fig. as well as 
for Fig. (denoted with the ^). All fluxes are in 
units of cm“^ s“^, and the column depths are in 
units of cm“^. 
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5.2.2. Scenarios involving surface sinks for O 2 


O2 has several potential sinks on an abiotic 
world, which we can cast in the form of mass bal¬ 
ance equations, following Kharecha et al. (20051. 
(Our full set of equations is listed in the Ap¬ 
pendix.) It is possible that high-temperature 
supercritical water (for example, as part of hy¬ 
drothermal circulation) could provide a medium 
in which CO and O2 can recombine. Hydrother¬ 
mal circulation through the midocean ridges would 
then act as a sink for both species through their 


direct reaction with one another (Helling & Tester 


19881. Oxygenated waters circulating through the 


midocean ridges could also oxidize new oceanic 
crust, creating an additional sink for O2. The 
timescale for circulation of ocean water through 
the midocean ridges on the modern Earth is ~10 
My (Isley 19951, and could have been lower on 
the early Earth if the geothermal heat flux was 
higher. If this is the only significant surface sink 
for O2, the O2 flux into the surface ocean must 
equal the flux through the midocean ridges. We 
calculate a lower limit of ~7xl0“® cm s“^ for the 
O2 deposition velocity in this case. If ydep{02) is 
set equal to this value, and if Vdep(CO)=10“® cm 
s“^, the calculated O2 mixing ratio at the surface 
is ~2xl0“^. The corresponding downward flux of 
O2 is ~10^^ mol yr“^. So, the midocean ridge sink 
for O2 would not prevent the occurrence of a false 
positive for life. 


A second potential surface sink for O2 is the ox¬ 
idation of formate in solution. If we assume that 
formate and O2 react in solution but that CO and 
O2 do not, and if this reaction is the only sink 
for both species, we can calculate the deposition 
velocity for O2 as a function of the partial pres¬ 
sure of CO (because CO controls the amount of 
dissolved formate). For a case with pco ~ 1 bar, 
Vdep(O2)~10“^^ cm s“^ - well below the deposi¬ 
tion velocity estimated for flow through the mi¬ 
docean ridges. The partial pressure of CO must 
be many orders of magnitude higher than this for 
formate to be a large sink for O2. That is unlikely, 
because it would require more carbon than is avail¬ 
able in the Earth’s surficial inventory. O2, how¬ 
ever, can be a substantial sink for formate. This 
would keep dissolved formate at relatively low con¬ 
centrations, thereby preventing formate from de¬ 
composing back into CO. In this case, an assumed 
CO deposition velocity of I0“® cm s“^ is realistic. 


Finally, it is conceivable that CO and O2 could 
react directly in solution, even though the rate 
constant for this reaction in a terrestrial ocean 
(where metal cations may catalyze the recombi¬ 
nation of CO and O2) is unknown. If this reaction 
were fast, dissolved CO and O2 concentrations 
should both approach zero, and their deposition 
velocities should both be close to their maximum, 
piston-velocity-limited value (more details in the 
Appendix). Fig. shows what the atmosphere 
would look like in the case for the planet orbiting 
GJ 876. The corresponding redox budget is shown 
in Table We have included a spectrum of this 
case in Fig. The surface O2 mixing ratio is 
^6x10“^, which is a factor of 4 below the Pro¬ 
terozoic O2 concentration inferred by |Planavsky| 


et al. (2014). The rapid recombination of dissolved 


CO and O2 would be enough to prevent an M-star 
false positive for life. 

In all of the cases examined so far, we have ig¬ 
nored another potentially important sink for O2, 
namely, reaction with dissolved ferrous iron in the 
ocean. On early Earth, the deposition of banded 
iron-formations suggests that the deep oceans con¬ 
tained abundant ferrous iron ( [Holland 1973). If 
this were true on the M-star planets simulated 
above, some of the O2 entering the ocean should 
react with Fe^’*' rather than with CO, thereby low¬ 
ering the concentration of dissolved O2 and in¬ 
creasing its deposition velocity. If the flux of fer¬ 
rous iron entering the oceans was sufficiently large, 
the dissolved O2 concentration should be zero, 
and the O2 deposition velocity should approach 
its maximum value of 1.4x10“^ cm s“^ (but the 
CO deposition velocity would not, so this differs 
from the case just discussed). Holland (2006) esti¬ 
mated the hydrothermal ferrous iron flux enter¬ 
ing the present oceans as 3x10^^ mol yr“^, or 
1.1x10^° Fe^"*" cm“^ s“^. Assume for now that 
the iron flux into the early oceans was the same. If 
the iron is oxidized to magnetite, as in eq. then 
l/6th this amount of O2 is consumed: 6 FeO -I- O2 
—>■ 2 Fe304. The maximum O2 flux that could be 
taken up by reaction with ferrous iron would be 
~2xI0® cm“^ s“^, or about 1% of the calculated 
downward flux of O2 in the I0W-O2 simulation. So, 
if Holland’s numbers apply to exoplanets, ferrous 
iron oxidation would be only a minor sink for O2. 

That said, ferrous iron fluxes on the early Earth 
could have been one to two orders of magnitude 
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higher than Holland’s estimate, depending on a 


number of pa rameters (Kump & Holland 1992 
Kasting|[2013 ). If all of the Fe^+ reacted with O2, 


and if the remainder of the O2 recombined with 
CO, then the dissolved O2 concentration should 
still approach zero, and CO would still be present 
in excess, keeping its deposition velocity close to 
the abiotic value of 10“® cm s“^. If CO is de¬ 
posited at this rate, and if O2 is deposited at its 
maximum rate, the atmosphere in our model en¬ 
ters CO runaway, in which the calculated CO con¬ 
centration continues to increase indefinitely. In 
reality, the CO concentration would eventually be 
limited by the finite source of CO2, and so the 
runaway would stop at some point; however, a 
dense, CO-dominated atmosphere might still build 
up. We have not explored such scenarios explic¬ 
itly; however, we should point out that this type 
of atmosphere could potentially be identified from 
the presence of CO absorption features centered 
at 2.35 and 4.7 /xm. These measurements would 
require a space-based telescope with near-IR capa¬ 
bilities extending to at least ^2.5 /im to capture 
the closest CO feature. The 4.7 fim CO feature 
would be more difficult to observe in direct imag¬ 
ing, since it overlaps with the 4.3 /xm CO2 band, 
and it occurs near the minimum of the combined 
planet-and-star emission. However, both the 2.35 
/xm and 4.7 /xm CO bands could be observable 
with transit transmission spectroscopy. 

Just prior to the onset of CO runaway, the at¬ 
mosphere just discussed closely resembles the low- 
02 atmosphere shown in Fig. At higher CO 
concentrations, calculated O2 levels fall. Thus, 
for all these types of atmospheres, O2 is below 
the Proterozoic O2 estimate of |Planavsky et al.| 

(pMIl). 


5.3. 


Comparisons between planets orbit¬ 
ing F, G, K, and M stars 


Tian et al. (2014) suggested that the high stel¬ 


lar far-ultraviolet to near-ultraviolet (FUV/NUV) 
ratio was responsible for the buildup of O2 in 
their M-star case, while [Domagal-Goldman et al.| 
(2014) pointed instead to the total FUV flux as 


the mechanism for building up O3 in their F-star 
case, as discussed in Section |5.2[ In Fig. [Tj we 
have explored this issue by arbitrarily changing 
the FUV/NUV flux ratio in the solar spectrum. 
We did this by decreasing the NUV flux and hold¬ 


ing the FUV flux constant (NUV is decreasing to 
the right in Fig. [^. The break between the FUV 
and the NUV is taken to be 170 nm. We have 
also plotted the FUV/NUV ratios for the other 
stars in our study. This was done while using our 
‘worst-case’ scenario (5% CO2, no rainout of CO 
and O2, Vdep(CO)=10“® cm s“\ Vdep{O 2 )= 0 . cm 
s“^). As the figure shows, the atmosphere transi¬ 
tions from low O2 to high O2 as the FUV/NUV 
ratio increases. Note that e Eridani (a K star) 
and a Bootis (an F star) have nearly the same 
FUV/NUV ratio, and that this ratio roughly cor¬ 
responds to the point of transition from low O2 to 
high O2. The K-star planet has substantially more 
O2 in its atmosphere than does the F-star planet, 
suggesting that the total UV flux is also an impor¬ 
tant factor. This effect appears to go the opposite 
way from that predicted by Domagal-Goldman et 
ah, as the K star has substantially lower FUV 
fluxes than the F star. We increased the FUV flux 
while holding NUV constant, to test this hypoth¬ 
esis (results not shown), and found that higher 
FUV increased O2 concentrations only slightly. 


The level of FUV has some impact on this 
result, as it is responsible for dissociating CO2, 
but the buildup of O2 is controlled principally 
by the NUV flux. We have not considered the 
effect of flare events, which has previously been 


examined for 02-dominated atmospheres (Segura 


et al. 2010), but such flares would likely only af¬ 


fect the atmosphere on short timescales, rather 
than the long-timescale steady state of our results. 
In steady state, Tian et al. suggested that the 
FUV/NUV ratio controls HO2 and H2O2 photol¬ 
ysis, which they suggest are sources of OH. We 
would point out, however, that these two species 
are actually just reservoirs for OH and that the 
source of OH is ultimately water vapor photoly¬ 
sis. Decreased water vapor photolysis, alongside 
the longer lifetime (and greater size) of the HO2 
and H2O2 reservoirs, limits the availability of OH 
for the recombination of CO and O2 in these at¬ 
mospheres. In our models, water vapor photolysis 
slows by a factor of 400 in the GJ 876 case com¬ 
pared to the solar case, and the size of the HO2 
and H2O2 reservoirs are 2-3 orders of magnitude 
larger. These factors together control the O2 con¬ 
tent in our model atmospheres. 


The results of our calculations are summarized 
in Fig. which shows vertical O2 profiles for plan- 
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ets around different types of stars. Solid curves 
represent ‘worst-case’ scenarios, that is, upper lim¬ 
its on abiotic O 2 ; the dashed curve for the M- 
star planet is our calculated O 2 concentration for 
a planet with an ocean in which O 2 reacts rapidly 
with CO. Planets with additional surface O 2 sinks 
would have even lower O 2 concentrations, due to 
the buildup of CO. For F- and G-star planets, the 
upper limit on O 2 is always well below the biotic 


(Proterozoic Earth) limit from Planavsky et al. 


(20141, so such planets should not exhibit false 
positives for life, according to our calculations. 
The upper limit on O 2 for the K-star planet is 
only marginally below the Planavsky et al. limit, 
suggesting that caution should be exercised when 
interpreting spectra from these planets. Or, to say 
this another way, one needs to carefully consider 
what the surface sinks for O 2 might be on such a 
planet. For the M-star planet, the upper limit on 
O 2 (near 0.1 PAL) is almost 100 times greater than 
the Planavsky et al. limit, and even our lower O 2 
estimate is close to that value. Thus, seeing O 2 
by itself on a planet around an M star may tell us 
little about whether that planet is actually inhab¬ 
ited. 


6. Discussion 

Thus far, we have focused on whether a calcu¬ 
lated O 2 or O 3 concentration could be higher than 
those on the Proterozoic Earth. But a second im¬ 
portant issue is whether or not we could detect the 
spectral features of abiotically generated O 2 and 
O 3 . Fig. 1^ shows synthetic spectra for the abiotic 
atmospheres with our ‘worst-case’ boundary con¬ 
ditions, i.e., our highest concentrations of O 2 and 
O 3 . We find that abiotic ozone (O 3 ) could be de¬ 
tectable for all cases except those with the Sun and 
e Eridani as host stars. The UV O 3 feature (the 
Hartley bands) is, unsurprisingly, most prominent 
for the GJ 876 ‘worst-case’ scenario (~ 6 % O 2 ). 
The results for AD Leonis are similar and are not 
shown. Curiously, the M-star planet atmosphere 
with the lowest O 2 concentration has deeper O 3 
absorption in its reflectivity spectra than the F- 
star planet; however the F-star planet produces a 
stronger O 3 signature in transit. This discrepancy 
can be explained by the distribution of O 3 in both 
atmospheres. While the I 0 W-O 2 , M-star planet 
has an O 3 column depth an order of magnitude 
greater than the F-star planet (see Table [^, the 


O 3 mixing ratio is significantly larger at high alti¬ 
tudes on the F-star planet (see Figures and |^. 
Because transit transmission observations will be 
insensitive to the lowest layers in the atmosphere 
in the UV and short-wavelength visible range due 
to Rayleigh scattering, the high-altitude O 3 in the 
F-star planet’s atmosphere would be more easily 
detectable. In neither of these cases should the de¬ 
tectability of O 3 in the UV be considered evidence 
for life. 

Our synthetic direct imaging and transit trans¬ 
mission spectra show that the O 2 -A band would 
be potentially detectable for M-star planet atmo¬ 
spheres defined by our ‘worst-case’ boundary con¬ 
ditions. For these cases, the O 2 -A band strength 
would be comparable to that of the modern Earth, 
as the band strength does not increase linearly 
with O 2 concentration at high O 2 levels. 

7. Conclusions 

Abiotic O 2 concentrations above ^10“^ PAL 
(the estimated O 2 concentration for the Protero¬ 
zoic Earth) could constitute a false positive for life, 
even if they are below the detection threshold of 
current and planned space and ground-based tele¬ 
scopes. O 2 concentrations this high, or higher, 
are readily produced on planets that lose large 
amounts of water. This includes runaway green¬ 
house planets within the inner edge of the HZ 
and M-star planets that are devolatilized during 
their pre-main-sequence evolution, which should 
be recognizable by the lack of water vapor in their 
atmospheres, and planets whose atmospheres are 
deficient in nitrogen and carbon (and where H 2 O 
is the major constituent), which should lack N 2 - 
N 2 absorption features. Frozen planets outside 
the outer edge of the HZ can also build up 02 - 
rich atmospheres, but these planets should also 
have H 20 -poor atmospheres. Thus, if the planet’s 
atmospheric composition is constrained, none of 
these false positives are likely to lead to significant 
confusion for future UVOIR exoplanet characteri¬ 
zation missions. 

Planets more like the Earth, be., with liquid 
water present on their surfaces and orbiting F- and 
G-type stars, are also unlikely to build up 02 -rich 
atmospheres, given reasonable assumptions about 
volcanic outgassing rates and surface oxidation. 
Earth-like planets orbiting K- and (especially) M- 


19 




type stars are a potentially bigger problem if they 
somehow retain surface water after the pre-main- 
sequence phase of stellar evolution. Such plan¬ 
ets are more prone to accumulate O 2 because of 
the higher far-UV to near-UV stellar flux ratio, as 
compared to the Sun. Photochemically produced 
abiotic O 2 levels that exceed the Proterozoic O 2 
estimates (< 10“^ PAL), such as those predicted 
by Tian et al. ( 2014| ), could constitute a possible 
false positive for life. However, high O 2 concen¬ 
trations are precluded for planets with sufficient 
surface sinks for O 2 (e.g., dissolved ferrous iron 
in the oceans), or if CO and O 2 react rapidly in 
solution. Laboratory experiments are needed to 
better determine the aqueous reaction rate of these 
gases in the presence of (catalytic) metal ions that 
would likely be present in planetary oceans. In the 
meantime, one should not rule out O 2 by itself as a 
potential biosignature, if caution and context are 
used appropriately. 
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Wavelength (nm) 


Fig. 1.— The top-of-atmosphere stellar flux for a 
planet located at 1 AU equivalent distance (receiv¬ 
ing 1,360 W of stellar radiation). The planets 
discussed here are located at 1.3 AU equivalent, 
which corresponds to a ~40% reduction in flux. 
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CO2 + hv 

0 + 0 + M - 

— CO + 0 
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^ O2 + M 


Troposphere 



H2O + hv ^ H + OH 


H + O2 + M ^ HO2 + M 
CO + OH ^ CO2 + H 
HO2 + O ^ O2 + OH 


CO O2 CO + O —" CO2 



3F^^+ V 2 O 2 + 60H'—^ Fe304 + 3 H 2 O 


2CO + 02^ 2 CO 2 


Fig. 2.— Summarizing the reactions that control 
abiotic O 2 levels in terrestrial atmospheres and 
oceans. 
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Mixing Ratio 


Mixing Ratio 


Fig. 3.— Comparative species profiles for the 
high-C02 cases shown in Fig. 6 of Hu et al. (2012) 
(left panels) and from our work (right panels). 
A 1-bar, 90 percent CO 2 atmosphere is assumed. 
Our model balances the global redox budget; Hu 
et al. balance only the atmospheric redox budget. 
Solid curves correspond to an H 2 outgassing rate 
of 3x10^° cm“^ s“^; dashed curves are for one- 
tenth that rate, and the dotted curves correspond 
to H 2 S outgassing only. The upper limit on Pro¬ 
terozoic O 2 proposed by Planavsky et al. (2014) is 
shown as a green line segment and arrow. 
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Mixing Ratio 


Mixing Ratio 


Fig. 4.— Results from our model for comparison 
with Figure 3 from Tian et al. (2014). Left-hand 
panel is for the Sun; right-hand panel is for the M 
star, GJ 876. We used their lower boundary con¬ 
ditions, outgassing rates, and ferrous iron fluxes 
(see text). We have modified the distance to the 
planet orbiting GJ 876 to be consistent with other 
results within this paper. The Planavsky et al. 
upper limit on Proterozoic O 2 is shown, as in Fig. 
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Fig. 5.— Using physically motivated boundary 
conditions (no sinks for O 2 , only abiotic sinks for 
CO), this figure shows the worst-case (highest O 2 ) 
scenario for abiotic planets. We include a Bootis 
(an F star), AD Leo (an M star), and e Eridani 
(a K star) for completeness. Note that the planet 
around AD Leo is similar to the one around GJ 
876 except that it has a higher abundance of CO. 
Again, the upper limit on Proterozoic O 2 is shown. 
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Mixing Ratio 


Fig. 6.— Realistic boundary conditions Vdep(02) 
= 1.4x10“^ cm s“^, VdepiCO) = 1.2x10“"^ cm 
s“^, no ferrous iron flux, normal rainout except 
no rainout of CO or O 2 . Only the GJ 876 (M- 
star) planet is shown, as the change in boundary 
conditions has little impact on either the solar or 
e Eridani planets, and as the effect on the AD Leo 
planet is nearly identical to that shown here. The 
upper limit on Proterozoic O 2 is shown. 
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GJ 876 



FUV:NUV Ratio 

Fig. 7.— Diagram showing the effect of changes 
in UV flux on the predicted surface O 2 concen¬ 
tration for a planet orbiting the sun. Decreasing 
the NUV (175 - 320 nm) portion of the solar spec¬ 
trum while holding the FUV (<175 nm) portion 
constant causes O 2 to increase. The boxes repre¬ 
sent the unmodified FUViNUV ratio for each of 
the stars in this study. 
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02 Mixing Ratio 


Fig. 8.— Calculated O 2 profiles for abiotic plan¬ 
ets orbiting a Bootis (‘F’), the Sun (‘G’), e Eridani 
(‘K’), and GJ 876 (‘M (range)’). Solid curves rep¬ 
resent upper limits on O 2 for cases where volcanic 
outgassing rates and surface O 2 sinks are small. 
The dashed curve represents the ‘I0W-O2’ case for 
the M-star planet. The upper limit on Proterozoic 
O 2 from Planavsky et al. (2014) is again shown. 
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121—^^^^^^^^^—I—.—.—.—.—^—.—.—^^—I 
0.755 0.760 0.765 0.770 0.775 

Wavelength [^m] 


Fig. 9.— Spectra of abiotic (5% CO 2 ) atmo¬ 
spheres. Top left: UV (0.15-0.4 /rm) synthetic di¬ 
recting imaging reflectivity spectra for a planet or¬ 
biting a Bootis (blue), the Sun (black), e Eridani 
(orange), GJ 876 with ‘I 0 W-O 2 ’ boundary condi¬ 
tions (dashed red), and GJ 876 with ‘worst case’ 
boundary conditions (solid red). The absorption 
from the O 3 Hartley bands is evident. Top right: 
synthetic reflectivity spectrum in the wavelength 
region of the O 2 -A band (0.755-0.775 /rm). Bot¬ 
tom left: synthetic transit transmission spectra in 
the UV wavelength range in units of the absorbing 
radius of the atmosphere in kilometers. Bottom 
right: synthetic transit transmission spectra of the 
wavelength region containing the O 2 -A band. The 
spectra of AD Leonis are similar to those of the 
high -02 (‘worst-case’) GJ 876 scenario and not 
shown here. See section 4.3 for a description of 
the models used to generate these spectra. 
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A. Ocean chemistry and its effect on deposition velocities for O2 and CO 


As discussed in the text, reactions of O 2 and CO in solution can have a controlling effect on their deposition 
velocities, and this, in turn, can determine whether atmospheric O 2 concentrations become high enough to 
constitute a false positive for life. Here, we describe our model for ocean chemistry on an abiotic planet. 
We acknowledge at the start that several important reaction rate constants are essentially unknown. Thus, 
instead of providing definitive answers, we focus on what needs to be better constrained in order to calculate 
realistic O 2 and CO deposition velocities. 


We begin from the ocean model of Kharecha et al. (2005). That model included both a surface and a 
deep ocean, and it did not include O 2 . We have simplified this model by combining the surface and deep 
ocean into a single ocean box. This assumption is reasonable because the time scale for deep ocean overturn, 
about 1,000 years on the present Earth, is shorter than nearly all the chemical reaction time constants in 
the model. This conclusion could change, of course, as new reaction rates are measured, but for now this 
simplifying assumption seems justified. 


Our ocean model keeps track of the concentrations of three dissolved species: CO, O 2 , and HCOO“ 
(formate ion). The suggested chemical reactions involving these species are shown in Section 4.4 (eqs. 8p0 1. 
Given this reaction set, we can write a set of mass balance equations for dissolved species concentrations: 


Vj,iCO){acoPco - [CO ])7 + h[HCOO-]^j 

z 

khyd[OH~][CO]zj 
Vp{02)iao2P02 - [02])l 
fcf°[C'0][02]z7 + [HCOO-][02\zi 

khyd[OH-][CO]z-f = k2[HCOO-]^-f 


+kZ^[CO][02]z^ + <S>mor[CO] 


+ ^MOr[02] 

+k3[HCOO~]z"f + k4[HCOO~]z^ 


(Al) 

(A2) 


+k^f^^~[HCOO-][02]zj +^mor[HCOO-] (A3) 

These reactions represent conservation of dissolved O 2 , CO, and formate, respectively. The units for these 
fluxes are mol yr“^. The units conversion factor, C, allows us to express dissolved species concentrations 
in traditional units of mol L“^. We assume that all the reactions with unconstrained rates are first order. 
The known rate constants and other assumed variables (e.g., pH) in our calculations are listed in Table 
Circulation through the midocean ridges is assumed to be a potential sink for all three species. 

With a few assumptions, we can use our model to duplicate the abiotic deposition velocity for CO obtained 
by Kharecha et al. If CO and O 2 react slowly in solution (A:£p —)■ 0), but HCOO“ and O 2 react quickly 
—>■ 00 ), so that HCOO“ does not decompose to CO, and if ocean cycling through the midocean 
ridges is slow, eq. |A1| reduces to: 


Vp{CO){acoPco - [C'0])7 = khyd[OH ][CO]z^ 
Rearranging, we can solve for [CO]: 

[CO] = acoPco - . ■, , —' / —WvtT 

1 “h \zf Vp\CO) 

Substituting this into eq. we can rearrange terms and solve for V(iep(CO): 


'^air \ 


l + khyd[OH ]z/vp{CO) 


(A4) 


(A5) 


(A6) 


Here, we have assumed that pCO/nCO = l/nair, consistent with a 1-bar surface pressure, as well as 
assuming that the surface ocean resembles that of the modern Earth (pH«8, T=298 K). For an Earth-like 
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planet, Vrfep(CO)« 10“® cm s“^. This is the deposition velocity for CO if its fastest sink is hydration to 
formate, and if formate itself is relatively short-lived so that it does not decompose back into CO. If formate 
is long-lived, the CO deposition velocity would be controlled by acetate formation, and its calculated value 
would be ^4x10“^° cm s“^. If dissolved CO reacts directly with O 2 by eq. then the CO deposition 
velocity could be considerably larger, up to its piston-velocity-limited value of 1.2cm s“^. These 
calculated rates are sensitive to the assumed surface temperature, and are a factor of ^4 higher (lower) for a 
15 °C higher (lower) surface temperature. Looking to eqs. Al][A3 we can use this same process to calculate 
the deposition velocity for other potential sinks for O 2 and CO, which we will detail below. The results of 
these calculations have been included in Section 15.2.21 


Dissolved O 2 and CO have three potential sinks: 1) flow through the midocean ridges, during which O 2 
and CO could react with each other or with the rocks, 2) reaction between O 2 and formate in solution (in 
which CO is assumed to have hydrated to formate), and 3) direct reaction between O 2 and CO in solution 
to regenerate CO 2 . We will consider each of these in isolation, which is appropriate if one particular sink is 
dominant. 


A.l. O2 and the Midocean ridges 

If the only sink for O 2 is circulation through the midocean ridges then we can write: 

Vp{02)iao2P02 — [02])l = ^mor[02]- (A7) 

Setting these two terms equal to each other, solving for [O 2 ], and incorporating an expression for the 
downward flux of O 2 analogous to eq. [^in the main text, we can determine a lower limit of Vdep( 02 )~ 7 x 10 “® 
cm s“^. The sink for CO from this same process would result in Vdep(CO)«6x 10“® cm s“^, assuming that 
CO either reacts with O 2 at high temperatures or undergoes the water-gas shift reaction (eq. |^. The latter 
case provides a physical mechanism to support our method for balancing the redox budget with a return flux 
of H 2 , as the calculated deposition velocity from this process is comparable to the assumed abiotic deposition 
velocity for CO. 


A. 2 . O2 and formate 


If formate and O 2 react quickly in solution (he., » 1 L mol ^ s ^), and if we further assume 


that this reaction dominates both the O 2 and formate systems, we can set the terms from eqs. A1 and A2 
equal to one another: 


Vp{02){ao2P02 - [02]h = k, 


HCOO- 


[HC00-][02]z^. 


(A 8 ) 


To determine pco and poj simultaneously, the system would need to evolve forward in time. We can, 
however, solve for the deposition velocity as a function of pco, with the deposition velocity represented as 
a fraction of the maximum deposition velocity (thus reducing the system to one variable). For example, 
to achieve 90% of the maximum O 2 deposition velocity through reactions with formate, pco > 10® bars, 
assuming the ocean is fully saturated in CO ([CO] = acoPco)- H the ocean were undersaturated with 
respect to CO, or if the reaction rate between O 2 and formate were slower, pco would have to be even 
larger. For our worst-case scenario with pco ~ 0.01 bar, Vdep(O 2 )« 10 “^^ cm s“^. Increasing pco to 1 bar 
increases the calculated deposition velocity to ~ 10 “^^ cm s“^. 


A. 3 . O2 and CO 

We can calculate the necessary reaction rate for the consumption of all available O 2 via the direct oxidation 
of CO (eq. [^. We do this by neglecting all other O 2 sinks and equating the flux of O 2 into the ocean with 
the oxidation rate of CO: 

Vp{02){ao2P02 - [02])l = /c£P[C0][02]z7- (A9) 

We can use our I 0 W-O 2 scenario for GJ 876 to establish a minimum for [CO] based on the CO flux into 
the ocean and its lifetime with respect to circulation through the midocean ridges. If ^dep(02) is 90% of 
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its maximum, we can constrain >4x10^ L mol~^ s“^. If the actual value for this reaction is lower 
than this, or if the reaction is not first-order in either reactant, then Vciep( 02 ) would fall from its maximum; 
however, a rate even two orders of magnitude slower maintains Vdep( 02 ) at 10% of its maximum, suggesting 
that a non-zero reaction rate for the direct oxidation of CO would be sufficient to draw down surface O 2 
below the Planavsky et al. limit. 


32 



Table 6 

Model variables and constants 


Variable Value Units Definition 


Vp(CO) 

4.8x10-3 

cm s~^ 

Vpi02} 

6.0x10-3 

cm s~^ 

o^CO 

S.xlO--* 

mol bar~^ 

0 O 2 

l.xlO-3 

mol bar~^ 

Px 

- 

bar 

[^] 

- 

mol 

Zs 

IxlO'^ 

cm 

Z 

4x10® 

cm 

C 

6.02x1030 

L mol~^ cm~^ 

d 

0.00374 

L cm~^ s~^ mol~ 

7 

1.61x1033 

L s cm“^ yr~^ 

khyd 

g{-10570/T+25.6) 

L mol~^ s“^ 

k2 

6.4x10-® 

s-i 

kz 

2.7x10-® 

s-i 

^4 

8x10-34 

s“^ 

L.CO 

^ox 

* 

L mol“^ s“^ 

^MOR 

1.4x103’ 

L yr-i 

uHCOO~ 

'^ox 

* 

L mol“^ s“^ 

T 

273-298 

K 

pH 

8 

- 

Psurf 

1 

bar 


piston velocity (=D/z) for CO 
piston velocity (=D/z) for O 2 
Henry’s law coefficient for CO 
Henry’s law coefficient for O 2 
surface partial pressure of species X 
dissolved concentration of species X 
Surface ocean depth 
Global average ocean depth 
Units conversion factor 

Conversion from mol yr~^ to molecules cm“^ s“^ 
Conversion factor (= C/d) 

Hydration of CO to HCOO“ 
Photodegradation of HCOO“ to CO (surface ocean) 
Acetate formation from HCOO“ 

Thermal degradation of HCOO“ to CO 
Direct oxidation of CO (eq. 

Ocean flux through the midocean ridges (e.g., eq. 
Oxidation of HCOO~ (eq. |lo[ l 
Range in ocean temperatures (necessary for k^yd) 
Ocean pH (necessary for [OH“]j^ 10~® mol L“^) 
Earth surface pressure 


* Unconstrained rate constant 
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B. Reaction Table 
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Table 7 

Reactions and rates. 



0("D) + MO + M 1.8x10- 

0(^D) + O 2 O + O 2 3.3x10 




Table 7 —Continued 



8.85x10 





Table 7 —Continued 



+ M ^ H 2 CO + OH 





Table 7 —Continued 






Table 7 —Continued 



lO 

(N ^ 

^ ^ 2 7 2 

-11 

-12 

1 

o 

Oi 1 ' 1 

(M o 2 

(JU 1 1 

CD o O 

Pi ^ -H 

X 

H H X ^ X 


Oi 

^ ^ 

o ^ 

CO 

« (M CO ] 1—1 

CO t-h 


o o 

CN CN 


_L hC ^ 

+ o o 

lO , u 

ffi + ffi 

u ffi + 

t u ^ 


^ ffi t 

0+0 

+ ffi + 

® ^ 
ffi O ill 


ffi 

1 

o 

o 

o 




+ 

o 

z; 

+ 



ffi 

o 

o 

+ 

+ 

o 

o 

CN 

o 

o 

ffi 

CO 

ffi 

CO 

X 

o 

ffi 

+ 

CO 

+ 

+ 

M 

CN 

O 

(M 

o 

+ 

o 

CO 

+ 

o 

o 

CO 

X 

o 

t 

o 

T 

o 

w 

ffi 

o 

t 

s, 

t 

CN 

t 

CN 

t 

o 

z 

t 


o 

O 

X 

K 

X 

(N 

+ 

+ 

ffi 

o 

+ 

+ 

o 

Cl 

o 

o 

+ 

+ 

lO 

lO 

ffi 

(M 

lO 

X 

ca 

+ 

lO 

CO 

CO 

o 

CO 

ffi 

O 

o 

X 




d 



d 





^ 1 ko=3xl0-3i •(T/300)“3' 

1 koo=1.5xlO“^^ -(T/SOO) 











Table 7 —Continued 


— 









o 

o 

1—f 

o 




o 

o 

4451 

T 

150C 

T 






1 




lO 





X 

X 

X 




o 

CO 

II 

1 




o 


w 

1 


CO 

1 

,—1 


CO 

1 

CO 

CO 

00 

1 

Nl 


CO 

O 

9 

00 

C-i 

cs 

1 

1 

1 

1 

M 

a 

o 

r—1 

Cri 

Pi 

O 

T—I 

s 

Pi 

s 

o 

T—I 

o 

oi 


o 

T—I 

o 




CN 


o 

o 

T—I 

O 

o 

T—f 

, c>. 

X 

tn 


X 

O! 

U2 

O! 

X 

t- 

t- 


X 




1 

1 

X 

X 

X 

X 

0) 

00 

< 

<C 

q 

< 

< 

< 

q 

1 

o 

1 

o 

1 

o 

CO 

o 

CO 

1 

o 

CO 

(N 

1 

O 

t-H 

o 

l-H 


q 

q 

q 

q 

4^ 

CN 



1—f 




CO 


1—f 


00 

1 

o 

1 

o 

X 

X 



t-H 












X 

X 

X 


o 


o 


00 

t-H 

I-H 





l-H 









q 

l>- 

CD 



1—f 

X 

1—1 

X 

lO 

l-H 

X 

o 

X 

CN 

CO 

ci 

i 

o 








o 

§ 

o 

o 


o 

+ 

CO 

CO 

o 

+ 


+ 

+ 

CO 

o 

CO 

O 

O 

+ 

CO 

ffi 

Z; 

CN 

CN 

ffi 

-0- 

t 

ffi 

-4- 

0- 

0- 


^ o 

o o + 

® K o 


O W K O 

^ O + + 

+ + o o 

O O g2 CO 

CO £2 W ffi 


O O 
+ + 


O O 

+ + 

CN (M 

ffi ffi 

o o 


K + 

O O t 

'I' o ^ 


O K 

o W ^ 

ffi + o 

+ £ + 

K + pH 

^ O ^ 

tot 
O t W 


CO CO ^ 

ffi ffi o 

o o 





Reaction Rate [i 




a 



®lh 

in P 

1 1 

Vr 


18( 

T 





(M 



O 

1 

l| 

o 

OT 

00 

00 

r . 

1 

1 

O 

o 

a> 

1 

1 


CO 

0^ 


o 

1 

o 

T—H 


(N 


(N 

1 

O 

T—f 

1 

o 


m 

O 

X 

X 

1 

1 

Q 

1 

1 

X 

X 

<< 

< 

7 

LO 

o 

00 

o 

1—H 


o 

1—H 

o 

CO 

o 



o 

1-H 


c6 

X 

X 

Q 

X 

X 

1—I 




X 



CO 

CO 

CO 

CN 





to 



CN 

CO 

oi 




oolh c^lh 


o ® W 

o + + 

+ O o 

® O o 

ffi ffi ffi 

o + + 

t ffi £ 

« o o 

o t t 

+ ffi o 

o ° + 

u + lO 

n ffi " 

0^0 

o 


CN ^ CO 

O ffi ffi 

- o o 


I CN O 

+ a 

W + ffi 

+ M _l_ 

ffi o o 


(N 

o 

CO 







o 


o 

o 


o 

ffi 

CO 

a:i 

o 

CO 

+ 

O 

o 

m 

o 



o 

+ 

CO 

w 

+ 

o 

iM 

ffi 

+ 

CN 

+ 

CO 

+ 

W 

+ 

CO 

oa 

+ 

CO 

CN 

+ 

CO 

(M 

CO 

+ 

CN 

+ 

CO 

CN 

+ 

o 

o 


ffi 

O 

ffi 

O 

ffi 

A. 

K 

ffi 

A. 


ffi 

CO 

X 


t t t 


K 

O 

ffi 

+ 

O 

+ 

O 

+ 

CN 

O 

o 

X 

+ 

CO 

CO 

CO 

+ 

+ 

CO 

oa 

ffi 

CN 

ffi 

(M 

ffi 

ffi 

CO 

ffi 

CO 

X 


O W o 

s + g 
+ s “ 


^ o 

+ o W o 

o + + + 

23 O O M 


t t t t t 

o" O o g 8 

+ z + + I 

(yj ^ cn CQ + 




Table 7 —Continued 






b- 

b- 



0 

• 

1 








1 

1 



T— I 

,.61 

1 








0 

0 



X 






C^l 



X 

X 

M 

0 

CO 

1 

05 

CM 




05 

1 

1 

0 

00 

t- 

1 

1 

+ 



1 

1 

1 

Pi 

1 

0 

T—f 

1 

0 


CD 

CO 

0 

1—1 

1 

0 

1 —I 

05 


0 

0 

1 —1 

1 

0 

T— ( 

c /3 

< 

X 

q 

X 

q 

U 2 

1 

0 

lO 

1 

0 

X 

0 

CD 

X 

10 



C 5 ^ 

(N 

1 

X 

CD 

X 

t- 

X 

05 


1—1 





1—I 

CN 

§ 


0 


I> 



CN 

K 

CN 

0 

0 

0 

0 

0 

+ 

+ 

+ 

+ 

CO 

X 

t 

0 

CO 

t 

0 

CO 

0 

CO 

t 

CN 

0 

1 

CN 

0 

T 

fO 

r5 

1 

CN 

0 

K 

w 


u 

+ 

+ 

+ 

+ 

CO 

CO 


CO 


+ + + 

22 M K 

ffi ffi g 


22 ^ -d 

+ ® + + + 


_ O O CO 

O CO CO ^ 

CO p 5 ffi 


+ + + O 

O CO CO + 


It! _L ^ j_ O ^ 

kS + 2;+0_l4- 

^ t ^ I ^ ffi 5? z; 

I S T S t t T t 

£ + gt°Oo° 

Offi,ffi + +^ + 

, ':7 . hH 


® 5 z 
z 


+ ® w 

z z £ 

t t I 
O i K 
+ + + 


z 2 ; K 
2 


N 2 H 3 + H-> NHz + NHz 2.7x10 






Table 7 —Continued 







Notes— References: 
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Sander et al.| ( |2011); ^ Baulch et al. 


[NIST]; ^ pYiedrichs et al.| (|2002| [NIST jT ^lPinto et al | (|1980l ) ; 


DeMore 


et al. 

Tfi 


(19991; 


(1995) [NIST 


Baulch et ai| (|1994|) [NIST 
■ ■ " 1 “ 


Ashfold et al.| (|1981|); ^|Tian et al. 
" fWagner fc Wardlaw] (|1988D ' ' 


leMore et al. 



(| 1992|) [NIST]; ^iTsang fc Hampson| (|1986|) 

^iQehlers et al. (2000|) [NIST]; ['[Yung fc 


^fPrasad fc Huntress Jr| 


1 


rouard et al. (1989); ^'Choi & Lin 


&: Herron J 1991|) [NIST ] ; ^fGannon et a l. J 2008 )^ 


(1997) 


[NIST 

arpl. 2" 


"Veyret fc Lesclaux (1981[) [NIST 


^iHumpfer 

( |2b05| ) [ NISTJ^ 


fZahnle (1986Huebner & Giguere 


et al. 


et al. 


Krasnoperov & Michael! ( 2004|) [NIST ; ^'Dillon et al.|(|2007 

) [NIST ; ^Sillesen et al. (1991“ NIST]; ^fLaidler 

& Wojciechowski|(|196l| NIST]; ^fDeMore et al. (1992); ^ 

Kastingl (|1990D; ^'Turco et al.|(1982|); ^fGraham 


1979^ ; ^ fYung fc DeMore ( 1982[); ^fHoinan fc Schweinfurth|(|1981[ ); ^ fGladstone| (|1983[) ; ^fBogdanchikov 
2004) [NIST 7^ Adachi et al. ( 1986[) ; ^'Harding et al. ( 2005[) [NIST], Wen et al. ( 1989|; ^' Beng et al. 


NIST 


] 1999[) [NIST] ; ^^tachnik fc Molina (1987) [NIST]; ^' Baulch et al.| (|1976|); ^ Atkinson et al. ( 2004 ) _ 

Gordon et al. I ( |1971[); ^^ark fc Lin| (|1999[) [NIST]; "^fFagerstrom et al.| ( |1995[) [NIST]; "^fGo Si fc Westberg l 
(119911) [NIST]; '^fStief & Payne] (|1976|); "^[Gehring et al.| (|1969|) [NIST]; "^iKuhn & Atreya] (|1979|); "^fKasting 
5 |Lenb et al.| ( |1972D ; ^ fWennberg et al.| ( |1994D ]NIST]. 
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“ These reactions are particularly sensitive to the chosen rate constant. 

^ These reaction rates take the form: k(M,T) = ko{T)[M]/[l+ko{T)[m]/ko 


.{T)]-0.6^[l+[logio[ko{T)[M]/k^{T)]f 


where ko{T) has units of cm® molecules ^ s ^ and koo(r) has units of cm® molecules ^ s 

° The photolysis rates presented here are taken from the uppermost layer in the model (^100 km) for 
the Sun (first value) and GJ 876 (second value). Caution: the rates in the upper atmosphere are not good 
indicators for rates in the lower atmosphere. 

^ These reaction rates take the form: k(M,T) = [ko{T)kooiT)[M]]/[ko{T)[M] + kaoiT)], where ko has units 
of cm® molecules”^ s“^ and koo(T) has units of cm® molecules”^ s“^. If FC is given, use instead: logio(k) 
= log io[koiT)/[l + ko{T)/k^{T)]] + {logio{FC)/[l + logio{ko{T)/k^{T)y], following [Wagner & Wardlaw 
( fl988l ). 

® Estimated rate constant or branching ratio. 

^ Modified products. 
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